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C.  A.  Angell 
Department  of  Chemistry 
Purdue  University 
Vest  Lafayette,  Indiana  47907 

Studies  of  the  supercooling  phenoaenon  in  water  are  hardly  a  new  activity, 
the  first  review  having  been  published  in  1775.^  However,  it  is  only  in  the 
past  decade  that  a  high  level  of  activity,  directed  particularly  at  the  physi¬ 
cal  properties  of  water  in  the  supercooled  regiae,  has  been  manifested.  Apart 
froa  the  obvious  aotivation  of  seeking  data  on  this  important  liquid  over  the 
maximua  range  of  conditions,  this  activity  has  been  stimulated  by  the  recogni¬ 
tion  that,  as  temperature  is  decreased,  certain  theraodynaaic  and  transport 
properties  appear  to  be  diverging  as  at  a  critical  point.  This  unexpected 
characteristic  has  naturally  led  to  an  intense  curiosity  as  to  the  nature  of 
the  cooperative  phenoaenon  involved  and  has  promoted  an  increasing  number  of 

experiaental  and  theoretical  studies  since  the  first  evidence  was  presented. 

(2) 

The  position  of  the  'field'  at  the  moment  is  an  interesting  one  with  many  of 
the  fundamental  properties  now  determined  with  acceptable  precision  (given  the 
difficulties  of  performing  measurements  under  conditions  of  metastability,  with 
the  requirement  of  small  sample  techniques  that  this  implies)  and  several  non- 
coapl iaentary  theoretical  models  in  circulation  to  explain  them. 

Since  this  is  the  first  review  of  this  subject  area  to  appear  in  ARPC,  we 
will  give  some  coverage  to  the  earlier,  and  more  scattered  studies  before 
developing  the  subject  systematically  to  its  current  position. 


THE  NORMAL  LIMIT  OF 


AND  THE  POSSIBILITY  OF  VITRIFICATION 


Although  the  principal  interest  in  the  supercooling  of  vster  lies  in  the 
extraordinary  temperature  dependences  manifested  in  its  physical  properties  as 

o 

temperature  decreases  belov  0  C,  these  measurements  can  only  be  performed  while 
the  sample  remains  free  from  nuclei  of  the  crystalline  phase  which,  once 
present,  will  necessarily  terminate  the  measurements  in  progress.  Hence  we 
must  first  discuss  those  factors  which  determine  where  in  temperature  a  given 
sample  will  convert  to  the  stable  phase. 

In  bulk  samples  the  mechanism  of  generation  of  the  first  element  of  the 
crystalline  phase  is  almost  always  heterogeneous  in  character,  i.e.,  the  ini¬ 
tial  crystal  commences  to  grow  on  the  surface  of  some  solid  impurity  present  in 
suspension  or  at  the  walls  of  the  containing  vessel.  The  heterogeneous  mechan¬ 
ism  short-circuits  the  needs  for  a  large  free  energy  fluctuation  sufficient  to 
overcome  the  difference  between  the  free  energy  per  mole  of  microcrystals  of 
the  order  10-20  £  in  diameter,  and  that  of  the  macroscopic  crystalline  phase 
which  is  in  equilibrium  with  bulk  water  at  0°C.  In  extremely  clean  water, 
nucleation  is  much  less  probable  hence  supercooling  more  extensive. 

Meteorologists  and  rainmakers  have  a  great  natural  interest  in  the  effi¬ 
ciency  of  various  impurities  in  causing  the  nucleation  of  supercooled  water 
droplets,  since  these  constitute  the  substance  of  high  level  clouds  and  their 
crystallization  is  often  the  essential  precursor  step  to  the  development  of 
raindrops.  Much  of  the  early  work  on  crystallization  of  supercooled  water  has 
therefore  come  from  meteorologists.  This  subjeot,  and  the  work  preceding  1962 
is  the  subject  of  detailed  exposition  by  Mason  in  his  book  Clouds.  Rain,  and 

Rainmaking .  ^  Important  contributions  were  made  by  Biggs,  ^  Mossop  "  and 

(3) 

Mason  et  al.  who  performed  systematic  measurements  on  water  samples  of 


different  volumes  contaminated  by  controlled  quantities  of  different  solid 


materials,  and  established  more  clearly  than  hitherto  a  logarithmic  relation¬ 
ship  between  the  volume  of  the  sample  under  study  and  the  temperature  to  which 
it  could  be  supercooled  before  crystallization  occurred.  The  relationship  has 
a  simple  statistical  basis  in  the  probability  that  the  sample  under  study  will 
contain  a  seed  of  sufficient  potency  to  allow  a  nucleus-producing  fluctuation 
to  occur  in  the  'standard'  time  of  observation.  This  follows  from  the  assump¬ 
tion  that  the  probability  of  crystallization  is  a  function  only  of  the  funda¬ 
mental  nucleation  rate  (not  valid  for  viscous  liquids)  and  that  the  nucleation 

rate  J  (nuclei  per  unit  volume  per  unit  time)  is  determined  by  events  which  are 

(3) 

random  in  time  and  space  leading  to  the  relations 

lnj  *  const  +  kAT 


(1) 


laV  “  const  +  kAT 

where  V  is  the  volume  in  which  a  nucleation  event  will  be  observed  in  unit  time 
at  undercooling  AT.  The  parameter  k  characterizes  the  shape  of  the  hetero¬ 
geneous  'seed  potency'  distribution,  while  the  constant  is  determined  by  the 
concentration  of  seeds  of  all  types,  which  may  be  characterized  as  the  'total 
dirtiness'  of  the  sample. 

Nucleating  impurities  seem  to  be  ubiquitous  in  chemistry  and  meteorology 
laboratories,  and  they  can  only  be  removed  by  very  laborious  and  painful  proce¬ 
dures.  Such  a  procedure  has  been  described  by  Wylie ^  and  permitted  him  to 
supercool  a  sample  of  approximately  1  cc,  i.e.,  a  truly  bulk  sample,  to  -31°C. 
This  feat  has  never  been  repeated.  Untreated  bul'.  (>1  cc)  samples  of  water 
subject  to  slow  cooling  in  the  absence  of  agitation  cannot  usually  be  cooled 
below  -12°C.  A  less  involved  procedure  for  cooling  smaller  samples  contained 
in  0.5  mm  diameter  capillaries  reproducibly  to  -34°C  has  been  described  by 


(5) 


Mo  s sop 


and  utilized  by  others  since.  The  essential  step  is  the  prolonged 


distillation  of  the  primary  water  sample  to  remove,  in  the  first  place,  nucle¬ 
ating  impurities  contained  in  the  laboratory  atmosphere  over  the  sample,  which 
are  swept  out  by  the  evolving  vapor.  It  is  also  essential  to  avoid  mechanical 
entrapment  of  droplets  of  the  boiling  liquid  in  the  issuing  vapor. 

Nucleation  also  frequently  occurs,  even  in  very  clean  water,  at  centers 
located  in  the  glass  containing  vessels. The  precise  nature  of  these 
crystallization  centers  is  not  known,  but  they  appear  to  occur  randomly  in  the 
surface  of  pyrez  and  silica  glass,  and  may  develop  over  a  period  of  time  due 
presumably  to  some  leaching  action.  Coating  a  glass  sample  container  with 
silicone  can  reduce  surface  nucleation  probabilities.^ 

The  density  of  crystallization-indncing  impurities  in  water  never  seems  to 

-12 

be  greater  than  one  effective  center  per  10  1,  so  that  samples  which  have 

been  subdivided  to  this  volume,  as  in  the  formation  of  5  pm  droplets  emulsions 
using  non-nucleating  surfactants, ^ 8  10^  seem  to  supercool  reliably  to  tempera¬ 
tures  essentially  the  same  as  those  found  for  similar  size  droplets  produced  in 
cloud  chambers  swept  free  of  all  heterogeneous  nuclei  in  preliminary  cleansing 
operations.  Such  samples  thus  seem  to  be  suitable  for  performing  meas¬ 

urements  on  water  to  temperatures  approaching  the  fundamental  limit  imposed  by 
the  'homogeneous'  nucleation  phenomenon.  *  '  In  the  latter,  microcrystals 

of  dimensions  sufficient  for  thermodynamic  stability  are  generated  by  the  fun¬ 
damental  entropy  fluctuations  characteristic  of  the  supercooled  liquid  itself, 
which  become  adequate  only  at  and  below  the  'homogeneous  nucleation  tempera¬ 
ture.'  The  homogeneous  nucleation  temperature  is  itself  a  function  of  sample 
size,  though  the  volume  dependence,  according  to  standard  nucleation 
theory,  '  '  is  much  smaller  than  for  heterogeneous  nucleation.  The 

relation  between  the  two,  as  best  it  is  known  at  this  time,  is  illustrated  in 

Figure  1  taken  from  the  book  of  Mason.  ^  Some  additional  recent  data  are  to 

(14) 

be  seen  in  an  article  summarizing  current  Russian  work  by  Skripov. 


The  lowest  tempera¬ 
ture  to  which  even  the 
smallest  and  purest  drop¬ 
lets  of  water  (produced 
in  cloud  chamber  experi¬ 
ments)  have  ever  been 
supercooled  seems  to  be 
-42°C,<2'5,11)  though 
scattered  and  controver¬ 
sial  reports  of  cooling 
to  lower  temperatures 
exist  in  the  literature. 

The  probable  reason  fox 
this  limit  will  be  high¬ 
lighted  in  the  discussion  *  (after  Massu,  ref.  3)  reproduced  by  permission). 

of  physical  properties  of  the  supercooled  liquid  to  be  given  below. 

All  of  the  above  measurements  pertain  to  observations  made  at  1  atm.  pres¬ 
sure.  Just  as  the  freezing  point  of  water  is  lowered  by  the  application  of 
hydrostatic  pressure,  so  is  the  temperature  at  which  the  nucleation  probability 
becomes  high.  Thus  a  question  of  interest  is  the  extent  to  which  the  nuclea¬ 
tion  temperature  of  water  can  be  depressed  by  exercise  of  the  pressure  vari¬ 
able.  This  question  has  been  answered  in  recent  years  by  two  groups  studying 
the  behavior  of  emulsified  samples.  Xans  and  Barnaud^1^  at  the  Universite 
des  Bas  Pyrenees,  studied  water  droplets  immersed  in  a  hydrocarbon  dispersant 

phase  using  a  lanolin  surfactant,  to  measure  the  nucleatiou  temperatures  to 

(17) 

pressures  of  80  MPa  (800  bar),  while  Eanno,  Speedy  and  Angell  using  a 

( 9) 

sorbitan  tristearate  emulsion  recommended  by  Rasmussen  and  MaeCenzie,  ex¬ 


Figure  1.  Homogeneous  and  heterogeneous  nucleation  . 
temperatures  for  water  as  a  function  of  sample  volume 


tended  the  measurements  to  0.3  GPa.  The  latter  measurements,  seen  in  Figure  2, 
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show  that  s  lower  limit  of  -92°C  is  retched  at  0.2  GPa  where  a  reversal  in  the 
pressure  dependence  of  the  nncleation  temperature  occurs  due  to  the  preferen¬ 
tial  nucleation  of  a  denser  form  of  ice,  ice(III). 

The  question  arises  whether 
it  is  possible  to  cool  water  so 
quickly  that  the  nucleation  event 
will  be  bypassed  and  a  vitreous 
state  of  water  obtained  as  a  re-  — 
suit.  Attempts  to  perform  this  £ 

feat  date  back  over  many  decades,  §. 

E 
•» 

but  except  for  one  unpublished  *~ 
report  on  the  apparent  amorphous 
condition  of  a  sample  of  water 
which  had  been  splat-quenched  onto 
a  sapphire  block, (D.  R.  China nn, 

0  1000  2000  3000 

private  communication,  quoted  in  Pressure  (bars) 

ref.  18)  no  succese  had  been  Figure  2.  Metastable  range  for  liquid  water 

in  the  region  0.01-300  MPa  defined  by  emul- 
claimed  natal  the  recent  work  of  sion  nucleation  temperatures.  (After  Kanno 

(19-21)  et  al.  ref.  17,  reproduced  by  permission). 

Brugeller  and  Mayer.  The 

latter  workers  reported  that  a  sample  of  emulsified  water  contained  an  amor¬ 
phous  fraction  when  jet-quenched  at  high  pressures  into  a  liquid  cryomedium 

such  as  liquid  propane.  The  evidence  presented,  which  seems  fairly  convincing 

(22) 

although  is  not  universally  accepted,  was  the  presence  of  an  exotherm 

during  reheating  of  the  quenched  material  which  occurred  at  160  K,  the  same 

(23) 

temperature  at  which  vapor-deposited  amorphous  solid  water  AS W  (which  is 

clearly  a  non-  crystalline  ms  irial)  cr  tallizes  during  steady  reheating. 


Temperature  (*K) 


(20  21) 

More  recently  these  sane  workers,  *  using  an  improved  technique,  have 
reported  the  vitrification  of  pure  water  jetted  through  a  10  micron  aperture 
under  pressures  of  40  mPa  into  various  cryomedia.  The  exact  state  of  this 
apparently  vitreous  material  is  a  little  ambiguous  at  the  moment  because  of 
uncertainty  concerning  the  role  that  the  cryomedium  itself  plays  in  the  struc¬ 
ture  and  properties  of  the  resulting  amorphous  phase.  The  reheating  behavior, 

for  instance,  shows  different  characteristics  depending  on  the  extent  to  which 

(20  21) 

the  cryomedium  has  been  removed,  e.g.,  by  pumping.  '  This  is  also  true  of 
the  X-ray  structure  of  the  material  which  shows  interesting  differences  from 
that  of  vapor—  deposited  ASW  which  evidently  cannot  be  attributed  entirely  to 
the  presence  of  elements  of  the  cryomedium. 

(24) 

In  independent  work,  Oubochet  and  McDowell  give  evidence  for  the  vit¬ 
rification  of  tiny  drops  of  pure  water  caught  on  hydrophilic  carbon  film  elec¬ 
tron  microscope  screens  and  plunged  into  a  liquid  cryogen,  e.g.,  propane  or 
ethane.  The  evidence  presented  is  the  amorphous  state  electron  dif fractogram 
<. ,.ained  from  the  product  and  the  crystallization  of  this  material  to  Ice  I 

c 

observed  at  —135-140  E  during  reheating.  A  critical  cooling  rate  is  observed 

in  their  work,  since  only  Ice  1^  is  found  if  the  film  falls  into  boiling  liquid 

nitrogen,  while  a  nitrogen  slush  yields  Ice  I  .  The  recrystallization  tempera- 

c 

ture  depends  on  time  in  a  manner  similar  to  that  observed  for  the  amorphous 
vapor  deposit  by  Dowell  and  Rinfret,  though  seems  to  be  lower  than  expected 
from  DSC  observations  carried  out  at  various  heating  rates. ^5,31)  por  in_ 
stance,  Dubochet  et  ji^4(b))  £^a(j  crystallization  occurs  after  100  sec  at 

140  E,  whereas  even  at  a  constant  heating  rate  as  low  as  1  deg  per  100  sec, 
crystallization  in  Skripov's  deposits did  not  occur  until  158  E.  There 
is  an  unresolved  problem  in  determining  the  extent  to  which  a  given  kinetically 
unstable  glass  is  nucleated  during  the  cooling  process,  and  when  the  density  of 


nuclei  in  the  material  is  large  enough  that  the  material  should  be  regarded  as 

microcrystalline  even  though  amorphous  to  X-rays  or  electron  diffraction. 

Heavily  nucleated  glasses  should  crystallize  more  readily  during  heating  than 

their  nucleus-free  equivalents.  Detailed  comparisons  of  the  quenched  droplet 

and  vapor-deposited  ASW  dif fractograms  are  awaited  with  interest;  partial 

incorporation  of  the  cryogen  in  clathrate  structures  could  help  explain  an 

apparent  reduction  in  kinetic  stability  compared  with  the  vapor  deposit  which 

(23) 

can  itself  be  studied  at  leisure  at  140  K. 

The  importance  of  the  liquid-vitrification  work  is  two-fold.  Firstly  and 

less  significantly,  the  retention  of  an  amorphous  material  on  quenching  an 

initially  liquid  state  implies  that  at  least  transiently  the  liquid  existed, 

though  not  necessarily  in  a  state  of  internal  equil librium,  throughout  the 

extended  temperature  range  between  -40°C  where  it  would  normally  crystallize 

and  the  temperature  of  final  vitrification  which  is  not  exactly  known  for  the 

cooling  schedule  employed.  Of  greater  importance  is  the  opportunity  which  now 

exists  to  resolve  the  following  important  problem  concerning  the  relation 

between  liquid  water  and  the  much  studied  vapor-deposited  ASW. 

The  latter  material,  which  is  obtained  by  deposition  of  water  molecules 

from  the  vapor  onto  a  cold  substrate,  and  which  has  been  much  studied  in  recent 
(23 ) 

years,  (but  which  falls  somewhat  outside  of  the  scope  of  this  review  since 

it  is  palpably  a  solid  rather  than  a  supercooled  liquid  state  of  water)  is 

(28  29) 

somewhat  mysterious  at  this  time.  Two  different  groups  ’  have  reported 
the  existence  of  a  glass  transition  temperature  in  the  range  130-140  K,  the 
apparent  increase  of  heat  capacity  at  which  leads  to  an  incompatibility  with 
the  known  excess  entropy  of  liquid  water  at  the  lowest  temperature  to  which  it 
has  been  studied.  Johari^*^  pointed  out  that  it  was  not  possible  that  a 
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liquid  state  with  the  observed  excess  of  heat  capacity  over  that  of  the  amor¬ 
phous  solid  could  arrive  at  0°C,  on  heating,  with  as  little  total  entropy  as 
liquid  water  is  known  to  possess  at  that  temperature.  The  argument  requires 
either  that  amorphous  solid  water  is  thermodynamically  unconnected  with 
ordinary  liquid  water  or  that  there  is  some  error  in  the  glass  transition  and 
heat  capacity  measurements.  It  now  seems  probable  that  the  latter  is  the  case, 
for  reasons  given  immediately  below,  and  the  jet-quenched  vitreous  material  may 
well  provide  the  evidence  to  substantiate  this  claim. 

In  an  attempt  to  improve  the  quantitative  understanding  of  the  glass  tran¬ 
sition  in  the  vapor— deposited  vitreous  material  MacFarlane  in  this  labora- 
(31) 

tory  carefully  deposited  material  into  a  differential  scanning  calorimeter 

(23) 

pan  at  ~90  K  under  appropriate  slow  deposition,  high  vacuum,  conditions  and 
affected  a  cold  transfer  to  a  Perkin-Elmer  DSC-2  differential  calorimeter.  Re¬ 
peated  efforts  including  cycling,  and  annealing  at  130-140  K,  using  the  highest 
sensitivities  of  the  DCS-2  have  failed  to  reveal  any  glass  transition  phenome¬ 
non  in  the  temperature  range  expected,  while  the  deposit  always  crystallized 

sharply  at  138-160  K  during  the  final  heating  scan.  Thus  the  observations  of 
(29) 

Sugasaki  et  aJL  are  either  made  on  a  material  with  different  properties  from 
that  deposited  in  MacFarlane' s  experiments,  or  are  the  result  of  some  experi¬ 
mental  artifact. 

( 32) 

MacFarlane  and  Angel 1  have  concluded  from  the  DSC  studies  that  the  va¬ 
por  deposited  material,  which  has  been  shown  to  be  a  rather  defect-free  random 
network  structure,  has  no  glass  transition  in  the  range  up  to  160  K  and 

cannot  be  expected  to  exhibit  one  at  higher  temperatures  because,  as  with  most 
splat-cooled  metallic  glasses,  the  deposit  relaxes  more  rapidly  towards  the 
crystalline  state  than  towards  the  internally  equilibrated  amorphous  state 
(though  a  shoulder  which  appears  consistently  on  the  low  temperature  side  of 
the  crystallization  peak  could  prove  to  be  an  enthalpy  relaxation  effect,  i.e., 
a  glass  transition,  competing  with  the  exothermic  crystallization).  In  support 
of  this  argument,  we  note  that  ice  nucleates  and  grows  from  glass-forming  aque¬ 
ous  solutions  at  a  lower  temperature  than  ASW  crystallizes,  despite  the  much 
greater  glass-forming  tendency  of  the  former.  The  failure  of  extrapolations 


from  binary  solntion  data  to  correctly  predict  the  relatively  high  temperature 
of  this  'intrinsically  unobservable'  glass  transition  must  be  due  to  the  same 
factors  which  cause  the  glass  transition  temperature  of  silica  and  zinc  chlor- 
ide  to  be  rapidly  depressed  by  addition  of  second  components,  viz.,  the 
disruption  of  the  relaxation-resistant  network  structure  by  the  network- 
incompatible  second  components.  If  the  heat  of  crystallization  of  the 
jet-quenched  water  samples,  when  free  from  cryoliquid  contaminant,  can  be  shown 
to  be  the  same  as  that  of  the  vapor  deposited  vitreous  water  then  the  proposed 
reassignment  of  the  glass  transition  temperature  to  some  value  above  160  K 
(which  is  in  principle  undeterminable  at  'normal*  heating  rates  hence  is,  in  a 
real  sense,  fatuous  in  concept)  will  be  confirmed. 

Summarizing  this  section,  liquid  water  can  be  conveniently  studied  in  its 
metastable  states  down  to  approximately  -40°C  at  1  atm.  pressure  decreasing  to 
-92°C  at  2  kbars  of  pressure,  and  an  extension  of  this  range  for  ultra-fast 
measurements  (not  yet  designed)  is  perhaps  possible  in  principle  for  measure¬ 
ments  performed  on  microscopic  particles  during  cooling  processes  at  rates  of 
7  -1 

the  order  of  10  deg  sec  .  A  continuity  of  non-equilibrium  states  between 
liquid  water  and  amorphous  solid  water  has  not  yet  been  positively  established. 


though  it  appears  more  possible  now  than  previously. 


In  the  following  pages  we  will  consider  the  various  measurements  which 
have  been  made,  not  only  on  water,  but  also  on  various  aqueous  solutions  on 
which  the  focus  of  attention  has  been  the  solute  modification  of  pure  water 
properties  and  implications  for  the  nature  of  water  itself.  These  properties 
will  be  considered  under  the  separate  sections.  Structural  Characteristics, 
Thermodynamic  Properties,  and  Transport  and  Relaxation  Properties.  In  the 
description  of  these  properties  we  will  repeatedly  take  note  of  the  functional 
forms  for  temperature  and  pressure  dependences  of  physical  properties  exhibited 
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description  of  these  properties  we  will  repeatedly  take  note  of  the  functional 

forms  for  temperature  and  pressure  dependences  of  physical  properties  exhibited 

by  both  pure  water  and  its  solutions,  though  we  will  postpone  a  detailed 

discussion  of  these  forms  to  a  concluding  section  in  which  the  various  models, 

aided  by  computer  simulations,  will  be  discussed.  Inevitably,  some  overlapping 

of  the  interpretive  material  into  the  experimental  sections  will  be  necessary 

in  order  that  the  experiments  which  have  been  performed  appear  generally 

coherent  in  their  design  and  sequence.  Considerable  use  will  be  made  of  recent 

(2)  / 34) 

review  articles  by  the  author,  Lang  and  Ludemann.  For  a  popular  level 

account  of  the  anomalies  of  supercooled  and  superheated  water,  see  ref.  35. 

Sampling  techniques 

Since  as  the  previous  section  has  indicated,  it  is  extremely  difficult  to 
free  bulk  samples  of  water  of  all  nucleating  impurities,  the  pursuit  of  knowl¬ 
edge  of  the  physical  properties  of  water  substantially  below  0°C  has  for  the 
most  part  depended  on  the  use  of  one  or  other  'small  sample'  technique.  The 
choice  of  technique  has  depended  to  various  extents  on  the  patience  and  ingenu¬ 
ity  of  the  individual  investigators,  but  for  convenience  can  be  summarized 
under  three  headings,  as  has  been  done  in  more  detail  in  Ref.  2.  The  refer¬ 
ences  given  are  examples,  and  are  not  intended  to  be  exhaustive. 

THIN  FILMS  (TWO-DIMENSIONS  MACROSCOPIC)  The  use  of  thin  film  techniques  has 

proven  appropriate  for  certain  spectroscopic  studies,  and  crystallization 

o  ( 2 ) 

has  been  suppressed  at  least  to  -31  c.  Some  controversial  reports  of  much 

(37  38) 

deeper  supercooling  exist  in  the  literature.  ' 
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CAPILLARY  SAMPLES  (ONE-DIMENSION  MACROSCOPIC)  This  technique  has  been  used 

for  certain  PVT  determinations  where  the  length  of  a  thread  of  water  is  con¬ 
veniently  measured  in  glass  capillaries/3^  415  and  also  for  spin-echo  diffu- 
sion(42)  and  viscosity (43>  measurements,  and  ESR  studies/44'  Sample  survivals 
down  to  -33°C  have  been  reported  several  times.  In  1-2  pm  diameter  samples, 
supercooling  to  below  -35°C  is  reported/44'455  but  much  caution  against  in¬ 
terference  of  surface  effects  on  the  physical  property  under  study  must  be 
.  ,  (45) 

exercised . 

DROPLET  SAMPLES  (NO  DIMENSIONS  MACROSCOPIC)  The  use  of  single  droplets  of 
diameter  -1  mm  down  to  -33°C  has  been  reported^5  but  the  common  droplet 
technique  has  involved  tho  use  of  emulsion  methods,  in  which  various  surfac¬ 
tants  have  served/8  105  A  popular  and  generally  successful  choice  has  been 

developed  by  Rasmussen  and  Mackenzie  using  span  65  (sorbitan  tristearate)  with 

(8) 

a  light  hydrocarbon  dispersant  phase.  These  have  served  for  extensive 

PVT^85  spin  lattice  relaxation^485  and  even  spectroscopic^4^5  studies.  At  the 
limit  of  the  droplet  technique  we  find  the  possibility  of  studying  'nanoscopic* 
(as  opposed  to  'microscopic')  samples  in  newly  developed  micro-emulsion  systems 
in  which  the  only  -OH  groups  are  those  of  the  water  itself /^0»51) 

Structural  Characteristics 

In  the  last  few  years  both  scattering  and  resonance  spectroscopy  studies 
of  the  structural  characteristics  of  water  have  been  extended  to  subzero  tem¬ 
peratures  using  thin  film,  capillary,  and  droplet  sample  techniques,  in  some 
cases  reaching  to  the  homogeneous  nucleation  temperature.  Bundles  of  silica 
capillaries  were  utilized  by  Bosio  £t  a_l  to  study  the  elastic  scattering  of 
neutrons  and  X-rays  respectively  from  samples  down  to  -20°C, 


from  which 
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isochoric  temperature  differential  atrnctnre  factors  (which  are  less  subject  to 
correction  errors  than  the  structure  factors  themselves)  were  obtained,  to  the 
structure  factor  for  comparison  with  known  behavior  in  the  normal  temperature 
regime.  In  this  manner  they  were  able  to  observe  sharpening  of  features  re¬ 
lated  to  the  D-D,  and  D-0  ordering  in  the  liquid  and  see,  through  the  oppo¬ 
sitely  directed  shifts  of  first  and  second  S(Q)  peaks,  the  consequences  of 
increasing  the  rigid  and  tetrahedral  hydrogen  bonding  developing  with  decreas¬ 
ing  temperature.  This  work  extends  to  lower  temperatures  an  initial  isochoric 
study  by  Egelstaff  et  al^^  Emulsion  samples  have  also  been  briefly  studied 
but  background  difficulties  have  made  resolution  of  the  results  unsatisfactory 
so  far  (J.  Teixeira,  private  communication). 

(54) 

Capillary  samples  have  also  been  used  by  Bosio  et  al  to  perform  small 
angle  X-ray  scattering  studies  and  thereby  to  obtain  information  on  the  in¬ 
crease  in  mean-squared  density  fluctuation  which  compressibility  measurements 
had  earlier  shown  set  in  at  the  lower  tempertures.  The  small  angle  X-ray 
measurements  give  specific  information  on  the  correlation  length  of  the  fluc¬ 
tuations,  which  prove  to  be  of  order  8  X  at  -20°C.  Reduction  in  the  mean- 
squared  density  fluctuations  on  addition  of  ethanol  were  also  demonstrated  by 
this  technique.  This  is  of  special  interest  because  viscosity  and  conductivity 
measurements  suggest  structuring  in  the  solution  is  increased  (see  below). 
Quasielastic  neutron  scattering  studies  have  also  been  carried  down  to  -20°C, 
and  the  results  will  be  referred  to  in  discussion  of  transport  properties 
below. 

A  number  of  studies  of  the  vibrational  spectra  of  water  under  supercooled 
conditions  have  been  reported  in  the  last  few  years.  Raman  spectra  for  water, 
D^O  and  their  solutions,  have  been  studied  by  several  groups at  tempera¬ 
tures  ranging  down  to  -24°C,  while  infrared  spectra  in  both  the  intermolecul ar 
(58) 


far  IR  region 


and  the  intramolecular  overtone  regions  have  been  examined, 


in  the  latter  case  reaching  as  low  as  -35  C  at  1  atm.  ' 


and  -60  C  at 


20  MPa. 


D'Arrigo  el  al 


have  reported 


the  polarized  Raman  spectra  of  pore  z'Dj" 

I 

water  in  the  range  95  C  to  -24  C,  using 

i 

a  remarkable  source  of  pure  water  pro-  ^iSCj- 

3  i 

duced  in  phials  by  the  firm  Angellini  ~  j 

~  i 

SpA.  Despite  its  large  -10  cc  volume,  ^  p0i_ 

l  i 

> 

these  water  samples  may  be  supercooled  ~ 

A  * 

repeatedly  to  at  least  —  24°C  while  re- 

50- 

I 

melting  precisely  at  0  C.  D'Arrigo 

I 

et  al  showed  that  whereas  spectra  in  the 


range  of  0-95“C  repeatedly  yield  the  ^  3CX  ^  St5kes^,n 
much-discussed  isosbestic  point  at 

Figure  3.  Polarized  Raman  spectra  of 

3370  cm-1,  the  spectrum  at  -24°C  devi-  water  in  the  temperature  range  95  to 

-24°C.  (After  D'Arrigo  et  al  ref.  55 
ates  markedly  from  this  pattern  (see  reproduced  by  permission). 

Figure  3) .  These  authors  considered  the  distinction  to  lie  in  the  presence  of 

heterophase  fluctuations  involving  microscopic  incipient  true  ice  nuclei 7||  The 

latter  interpretation  is  not  directly  supported  by  the  more  detailed  study  of 

Bansil  e_t  al ^ ^  who  examined  both  the  3400  region  and  the  low  frequency  region 

o 

at  temperatures  down  to  -20  C,  though  increases  in  the  temperature  dependences 

of  the  'bonded'  OD  stretching  band  were  noted  in  the  supercooled  region.  An 

extrapolation  of  the  bandwidth  to  the  temperature  of  140  K,  the  expected  T  , 

8 

/ 53  50  50) 

where  data  on  the  vapor-deposited  amorphous  solid  are  available'  ’  * 

showed  close  concordance  (2450  cm  *  FTHH  *  80cm  *  vs.  2447  cm  *  FWHH  *  60 
cm  *,  respectively)  which  argues  in  favor  of  a  continuity  of  structural  states 
between  liquid  and  amorphous  solid  (though  according  to  MacFarlane's  and 


Angell's  recent  conclusions,  see  earlier  section,  the  comparison  should  be 
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made  at  a  higher  temperature.  Bansil  et  al  argued  for  two  classes  of  hydrogen 
bonds  (not  water  molecules),  strong  and  weak,  on  the  basis  of  their 
observations. 

A  subsequent  careful  study  of  the  composition  dependence  of  the  3400  cm  * 
region  for  H^O  in  dilute  solution  in  D^O  by  Wiaufe-Akenten  et  el^*^  has  iden¬ 
tified  the  low  frequency  component  of  both  H.O  and  D.O  Raman  spectra  as  a 

«  ** 

feature  determined  by  intermolecular  coupling.  The  weakening  of  the  coupling 
for  more  weakly  hydrogen  bonded  molecules  leads  to  the  high  frequency  feature 
(see  Figure  3) .  These  results  were  considered  more  consistent  with  a  continu¬ 
ous  distribution  of  hydrogen  bond  strengths  in  the  total  structure,  i.e.,  with 
a  continuum  model  for  the  liquid. 

Polarized  and  depolarized  Raman  spectra  for  the  same  frequency  region  have 
been  reported  by  Yeh  et  al,  who  also  studied  the  low  frequency  175  cm  1 
band  in  the  depolarized  spectrum.  This  band  is  believed  to  be  an 
intermolecular  H-bond  stretching  mode.  A  strong  increase  in  its  intensity 
observed  with  temperature  decrease  in  this  work  was  interpreted  in  terms  of 
water  molecule  cluster  (H^O)j  build-up  in  view  of  the  symmetry  requirement. 

Overtone  IR  studies  have  been  made  up  to  200  MPa  using  emulsion  samples  of 

both  pure  H^O,  and  H^O-D^O  mixtures,  only  the  former  giving  adequate  spectra  in 

(49) 

the  high  pressure  regime.  From  the  former,  Angell  and  Rodgers  were  able  to 
obtain  information  on  the  uncoupled  overtone  bonds  and  to  show  that  the  band 
cosuaonly  attributed  to  'weak'  or  'broken1  hydrogen  bonds  decreases  to  very 
small  intensities  at  the  low  temperature  extreme  leaving  behind  a  broad  struc¬ 
tured  band  which  could  be  attributed  to  the  fully  bonded  or  vitreous  structure 
and  had  features  in  common  with,  though  shifted  to  high  energies  from,  the 
corresponding  ice  1^  band.  This  band  compared  well  in  shape  with  the  spectrum 
of  i^O-glycerol  glasses  observed  at  -100°C.  The  results  could  not  be  inter¬ 
preted  within  the  framework  of  a  simple  two-state  model,  but  seemed  consistent 
with  the  Stillinger-Rahman  concept  of  'exciting  across  the  centroid,  *  the 


X—  .jl. 
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latter  being  a  static  population  of  H-bonds  of  intermediate  strength.  Increase 
of  pressure  was  found  to  reduce  the  enthalpy  change  associated  with  this  pro¬ 
cess.  (49(b))  gowever/  contrary  to  expectations  from  two-state  models,  the  num¬ 
ber  of  broken  bonds  does  not  increase  with  increasing  pressure  except  perhaps 
at  the  lowest  temperatures.  This  would  imply  that  the  strikingly  anomalous 
consequences  of  increasing  pressure  on  water  come  about  not  through  'breaking' 
hydrogen  bonds  but  through  redistributing  the  bonds  already  present  so  as  to 
have  fewer  groupings  of  strong  bonds  which  are  responsible  for  the  low  density 

fluctuations  (see  Interpretations  section  later). 

( 58) 

In  the  far  infrared,  L.  Boehm  has  shown  that  decrease  of  temperature 

to  -12°C  (or  to  -60°C  using  H^O  +  propylene  glycol  solutions),  results  in  a 

very  considerable  sharpening  and  shift  to  lower  frequencies  of  the  interesting 

feature  at  170°C  cm  1  reported  by  Williams  et  al/^  in  the  normal  temperature 

regime.  The  band  which  is  usually  assigned  to  a  translational  quasi-lattice 

mode,  is  also  observed  clearly  in  Raman  scattering  studies,  and  Bansil  et 
(63) 

al  have  recorded  temperature  effects  essentially  the  same  as  those  of 

Boehm. 

Thermodynamic  Properties 

PVT  MEASUREMENTS  The  outstanding  anomaly  characteristic  of  water  in  its 
stable  region  is  the  density  maximum  at  4°C,  and  it  has  been  a  natural  since 
earliest  times  to  enquire  into  the  behavior  of  the  density  in  the  negative 
expansivity  region  when  measurements  are  extended  below  0°C.  The  first  re¬ 
ported  measurements  date  back  to  1837  when  Despretz  performed  accurate  density 
studies  to  -9°C.*^  Using  capillary  samples,  such  measurement  .  have  been 
extended  more  recently  to  -34°C^^  while  emulsion  methods  have  pushed  the  low 
temperature  limit  down  to  -37°C.^  The  results  indicate  that  supercooling 
water  is  rapidly  approaching  the  low  density  of  ice  at  the  lowest  temperatures 
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of  measurements,  and  reasonable  extrapolations  indicate  that,  had  the  crystal¬ 
lization  event  not  occurred,  the  ice  and  water  density  plots  would  intersect  in 
the  vicinity  of  -(45-50°C),  see  Figure  4. 

The  peculiarity  of  this  behavior  is  best  displayed  through  the  derivative 
property,  the  thermal  expansivity,  which  is  found,  at  -34°C,  to  be  more  than 
three  times  as  negative  as  it  ever  becomes  positive  at  higher  temperatures  (ex¬ 
cept  near  the  critical  point).  The  behavior  of  the  density,  and  the  consis¬ 
tency  of  the  expansivity  measurements  in  the  supercooled  region  between  differ- 

(9) 

ent  capillary  and  emulsion  measurements,  and  their  overlap  with  those  of  the 
bulk  material,  is  well  demonstrated  in  Figure  4. 

The  expansivity  measure¬ 
ments  for  D^O  were  extended  to 

.  .  .  .  (41.65) 

high  pressures  by  Kanno, 

who  found  a  systematic  shift  of 
the  density  maximum  (at  11°C 
for  normal  pressure)  to  lower 
temperatures  with  increasing 
pressure.  At  the  limit  of  the 
study  the  TMD  occurred  at  -24°C 
at  120  MPa.  These  measurements 
have  recently  been  extended  to 
negative  pressures  of  the  order 
of  20  MPa  in  a  clever  experi¬ 
ment  by  Henderson  and  Speedy,  using  a  Bourdon-type  tube  technique.  At  -20 
MPa  the  density  maximum  for  B^O  occurs  at  7°C,  continuing  the  pattern  observed 
at  positive  pressures.  Extensions  of  these  measurements  using  remarkable 
microcapillary  spirals  have  permitted  observations  to  -18°C  under  negative 
pressures,  (Henderson,  B.  and  Speedy,  R.  J.,  private  communication). 


Figure  4.  Density  and  expansivity  dependence 
on  temperature  in  normal  and  supercooled  states 
of  H^O  and  D.,0  (from  ref.  2  by  permission). 


Comparable  anomalies  in  the  compressi¬ 


bility  were  recorded  by  Speedy  and 

Angell  at  low  pressures and  extended  to 

pressures  of  the  order  of  200  MPa  by  Eanno 

and  Angell. The  variations  follow 

closely  the  course  predicted  by  polynomial 

extrapolations  of  precise  measurements  in 

the  stable  region,  and  have  recently 

been  reconfirmed  by  new  extrapolations  of 

PVT  data  obtained  recently  using  a  novel 

piezo  thermal  method  by  Ter  Minassian  et 

alf67)  Their  Figure  5  is  reproduced  here.  Figure  5.  Compressibility  behavior 

of ‘water  up  to  200  MPa  by  extrapo- 

The  heat  of  compression  data  of  these  an-  lation  from  stable  range,  compared 

with  experimental  data  of  ref.  65. 

thors  extend  the  knowledge  of  water  to  600  (After  Ter  Minassian  et^  al,  ref.  67, 

reproduced  by  permission.) 

MPa  at  temperatures  ranging  down  to  -28  E. 

It  was  the  finding  that  the  low  pressure  compressibility  and  other  proper¬ 
ties  could  be  fitted  to  a  critical  type  equation 

X  -  Ax(T/T#-l)"T  (2) 

with  a  singular  temperature  just  below  the  homogeneous  nucleation  temperature 
which  first  drew  general  attention  to  the  interest  potential  of  water  in  the 
supercooled  region.  The  singular  temperature  obtained  by  constrained  fit¬ 

ting  of  a  critical  point  equation  to  the  compressibility  dataproves  to  follow 
closely  the  homogeneous  nucleation  temperature  over  a  range  of  pressures, 
essentially  as  shown  later  (in  Fig.  7)  from  the  analysis  of  NMR  relaxation  data. 

It  seems  rather  clear  from  these  relations  that  there  is  a  close  connection 
between  the  suggestive  temperature  dependences  of  these  intensive  thermodynamic 
properties  and  the  eventual  crystallization  event,  though  the  basis  for  the 
connection  remains  at  this  time  a  matter  of  some  controversy. 


The  adiabatic  compressibility  has  been  determined  by  ultrasonic  measure¬ 
ments  to  -26°C  by  Bacri  and  Rajonar ison^^  at  925  MHz  and  extended  to  -34°C  by 
(47) 

Trinh  and  Apfel,  using  an  ingenious  levitated  droplet  technique  and  54  KHz. 

waves.  A'_  determined  in  the  latter  work  shows  an  interesting  maximum  value  at 

about  -32°C,  a  feature  which  is  mimicked  at  slightly  higher  temperatures  by 
(69) 

C^.  Both  variations  make  a  striking  contrast  with  the  constant  pressure 

properties.  No  sign  of  the  maximum  is  seen  in  the  925  MHz  study.  Hyper¬ 
sonic  (55  GHz)  measurements  by  Teixeira  and  Leblond^^  using  the  Brillouin 
scattering  technique  down  to  -20°C  accord  better  with  the  925  MHz  velocities, 
but  show  a  different  absorption  a/f^  suggesting  a  dispersion  due  to  structural 
relaxation  may  be  at  hand.  Further  work  is  needed  in  this  area.  K  maxima 

like  that  observed  by  Trinh  and  Apfel  have  been  observed  in  H^O  +  ethanol 

(71) 

solutions  by  Conde  et  al .  The  ultrasonic  studies  have  been  the  subject  of 

a  detailed  review  by  D'Arrigo.^2* 

HEAT  CAPACITY  The  heat  capacity  of  supercooled  water  is  a  property  which 

lends  itself  to  measurement  by  the  emulsion  technique,  and  for  this  reason  its 

value  is  known  over  the  widest  temperature  range  of  all  the  intensive  thermo- 

(73) 

dynamic  properties.  Early  measurements  by  Rasmussen  et  al  and  Angell  et 

ai(18)  jjave  recently  been  refined,  and  correlated  with  similar  striking 

variations  encountered  in  the  superheated  regime — an  observation  of  which  more 

will  be  made  later  in  this  review.  Combination  of  C  with  the  above  PVT  data 

P 

has  yielded  the  constant  volume  heat  capacityC  which  has  a  maximum  at 
-23°C.<*9) 


OTHER  THERMODYNAMIC  PROPERTIES  The  vapor  pressure,  which  is  at  the  same  time 
the  most  fundamental,  and  the  longest-known  of  the  deeply  supercooled  water 


20 


properties  has  been  on  record  since  Regnault's  pioneering  work  to  -32.8°C  in 
(75) 

1847,  and  recently  confirmed  by  less  extensive  measnroments  of  Bottom- 

(76) 

ley.  Since  it  is  only  in  the  second  derivatives  of  the  free  energy  that 

the  anomalons  character  of  supercooled  water  becomes  clearly  displayed,  the 

vapor  pressure  measurements  themselves  show  little  to  excite  the  curiosity. 

The  remaining  thermodynamic  property  of  which  some  knowledge  is  available 

in  the  supercooled  regime,  is  the  low  frequency  relative  permittivity.  This 

(77) 

was  first  determined  using  an  emulsion  technique  by  Hasted  and  Shahidi  and 

quickly  supported  with  emulsion  measurements  using  a  different  analysis  by 
(78) 

Hodge  and  more  recently  by  the  co-axial  resonator  spectrometer  technique  of 

Bertolini  e_t  al  who  reached  -21°C  using  a  non-emulsion  sample  of  approximately 

(79) 

0.1  cc  volume.  Hermiz  (J.  B.  Hasted,  private  communication)  has  recently 

found  some  reason  to  doubt  the  viability  of  the  initial  emulsion  measurements, 
though  the  results  are  well  supported  by  independent  measurements  using  dif¬ 
ferent  techniques.  This  reviewer  is  unable  to  account  for  the  difficulties  in 

(79) 

emulsion  preparation  conscientiously  reported  by  Hermiz  in  this  connection. 

SOLUTION  MEASUREMENTS  It  is  not  possible  to  review  in  full  detail  the  find¬ 
ings  of  the  numerous  binary  solution  studies  which  have  now  been  reported  but 
we  will  summarize  briefly  their  intentions  and  results.  To  do  this  we  must 
first  introduce  some  of  the  basic  ideas  which  have  been  proposed  to  understand 
the  origin  of  the  curious  thermodynamic  characteristics  reviewed  above. 

Common  to  all  interpretations  is  the  notion  that,  with  decreasing  tempera¬ 
ture,  the  hydrogen  bonding  interactions  amongst  the  water  molecules  become  more 
effective  and  the  4-coordination  requirement  established  by  the  ice  rules 
results  in  an  increasing  open  network  character  (commonly,  though  in  some  ways 
misleadingly,  referred  to  as  'ice-likeness')  being  generated  in  the  molecular 


assemblage.  It  is  in  the  next  order  of  refinement,  namely,  in  specification  of 

the  types  of  configurations  encountered  at  the  microscopic  level,  that  the 

( 9  72) 

different  interpretations  show  their  distinctions.  One  class  of  model  ' 
considers  the  structured  groupings  to  be  characterized  by  actual  crystalline 
ice-like  arrangements  of  molecules  describable  as  Frenkel  heterophase  fluctua¬ 
tions,  and  generating  in  effect  ultramicrocrystals  from  the  otherwise  water¬ 
like  assemblages.  A  second  class insists  less  on  the  symmetry  of  ice 
but  expects  that  the  low  density  fluctuations  have  the  volumetric  charac¬ 
teristics  of  ice,  while  allowing  the  arrangements  to  be  spatially  random  in 
accord  with  the  'random  network'  concept.  A  third  group  empha¬ 

sizes  the  importance  of  5-membered  ring  groupings  such  as  are  exhibited  in 
clathrate  hydrate  compounds  but  not  in  the  ices.  Originally  conceptualized  by 
Pauling  as  a  favored  arrangement  in  stable  liquid  water,  this  idea  draws  some 
support  from  the  manner  in  which  such  groupings  can  introduce  a  cooperative 

characteristic  into  the  structuring  of  water  with  decreasing  tempera- 

t  (85.86) 
ture. 


With  these  basic  notions  now  introduced,  one  can  discuss  the  motivation 
for  solution  property  studies.  These  are  several-fold. 

An  initial  venture  has  been  based  on  the  idea  that,  since  the  open  struc¬ 
tures  are  determined  by  the  peculiar  charge  distribution  around  the  water 
molecule,  they  will  be  interrupted  by  the  presence  of  hydrophilic  solutes  such 
as  ions  and  other  strongly  hydrogen  bonding  molecules.  It  is  then  supposed 
that  such  binary  solutions  may,  at  sufficient  concentrations,  provide  a  basis 
for  deducing,  by  extrapolation  to  pure  water,  the  properties  of  water  in  the 

absence  of  framework-building  water-water  interactions.  This  has  been  the 

(69  87) 

particular  focus  of  Oguni  and  Angell  '  in  their  attempt  to  separate  an 


'anomalous'  component  from  the  total  thermodynamic  properties  in  order  to 
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W 


better  assess  the  appropriateness  of  critical  laws  for  the  'anomalous'  part  of 

the  thermodynamic  properties.  Startling  decreases  in  C  on  solnte  additions 

P 

( 69) 

are  observed  in  the  deep  supercooling  region. 

( 85 ) 

Stimulated  by  the  suggestion  of  Stillinger,  there  has  been  consider- 

( 88— 91) 

able  activity  in  seeking  solution  information  on  the  behavior  of  hydro- 

phobic  solutes  around  which  the  water  molecules  can  build  clathrate  cages. 

Such  induced  cage-building  could  be  viewed  as  a  positive  structuring  influence 

in  which  case  it  is  anticipated  that  the  singular  state  would  be  reached  more 

rapidly  with  decreasing  temperature  and  the  anomalous  behavior  therefore  en- 

( 92 ) 

hanced.  Indeed  early  studies  by  Vada  had  indicated  that  the  temperature  of 

the  density  maximum  of  water  was  initially  increased  by  incorporation  of 

ethanol  into  the  system.  These  ideas  have  been  tested  using  thermodynamic 

(88) 

measurements  by  Oguni  who  has  studied  the  variation  of  the  homogeneous 

nucleation  temperature  for  ice  with  solute  concentration  for  many  solutes,  and 
most  recently  by  Sorensen, who  studied  the  partial  molar  volume  of  solu¬ 
tions  with  hydrazine  and  t-butyl  alcohol  as  solutes  at  temperatures  down  to 
-30°C.  Except  for  Wada's  observation,  no  increases  in  the  anomalous  thermo¬ 
dynamic  characteristics  have  been  observed  (however,  see  Transport  Properties). 


Indeed  decreases  in  iC  as  abrupt  as  those  observed  for  C  with  all-hydrophilic 

a  p 

rved  for  EtOH  additions  by 

(7^  an  interesting  maximum  in  Ki  vs.  T  was  observed 


solutes^111^  have  been  observed  for  EtOH  additions  by  Conde  et  £l  (by  light 


scattering  measurements: 

in  solutions  of  >5%  EtOH  in  this  work).  However,  the  clathrating  concept  has 

( 88—91 ) 

survived  since  it  has  been  recognized  that  the  actual  mean-squared 

density  fluctuations,  which  determine  the  magnitude  of  kL,  a,  and  C  can  be 

T  P 

damped  by  including  guest  molecules  into  cages. 


It  is  in  the  transport  properties  of  supercooled  water  that  we  find  some 
of  the  most  unusual  behavior,  and  also  some  of  the  best  evidence  for  the 
reality  of  the  putative  singularity  underlying  the  behavior  of  this  extraordi¬ 
nary  liquid. 


VISCOSITY  The  most  familiar,  though  not  the  most  fundamental,  of  the  trans¬ 
port  properties  is  the  shear  viscosity,  which  at  1  atm.  pressure  has  been  now 
studied  over  a  temperature  range  extending  from  ISO  to  -3S°C  with  a  high  degree 
of  concordance  amongst  the  different  investigations.  The  earlier  measurements 
in  the  supercooled  state  made  by  Hallett^^a^  using  a  cold  capillary 
viscometry  technique,  have  recently  been  extended  to  -3S°C  by  Osipov  e_t  al^^ 
using  a  briefly  described  moving  slug  method.  The  accuracy  of  the  data  at  the 
lower  extreme  appear  to  be  about  ±1%  and  the  entire  temperature  range  can  be 
fitted,  (see  Figure  9  below)  to  the  singular  eqnation,  Eq.  (2)  above,  with  a 
singular  temperature  of  228  K  in  excellent  agreement  with  the  most  recent  and 
probably  reliable  determination  made  by  electrical  conductance  methods  to  be 
described  below. 

More  striking  than  the  highly  non-Arrhenius  temperatnre  dependence  (which 
after  all  is  seen  qualitatively  in  many  viscous  liquid  systems,  though  in 
detail  fitting  a  different  equation  form)  is  the  anomalous  presure  dependence 
of  the  viscosity.  Though  not  yet  studied  below  0°C,  the  existing  data  sup¬ 
ported  by  results  of  diffusivity  and  relaxation  time  measurements  at  lower 
temperatures,  show  clearly  that  a  striking  negative  presssure  dependence  of  the 
viscosity  develops  as  the  temperature  falls  below  20°C.  Such  an  anomaly  is 
actually  to  be  expected  from  theories  which  link  the  transport  properties. 


i.e.,  the  rate  of  exploration  of  the  system's  phase  space,  with  the  configura- 

(93) 

tional  entropy  of  the  liquid,  since  a  Maxwell  relation, 

l1)S\  .  _/hV\  fo 


(IS')  =  -(M\ 

V-&r'T  '■yr'r 


shows  S  increases  with  P  in  a  negative  expansivity  regime.  Studies  of  the 
viscosity  at  subzero  temperatures  and  pressures  to  200  MPa  are  currently  being 
performed  by  M.  Rogalski  in  the  author's  laboratory  using  a  micro-falling  ball 
(mercury  bead)  viscometer  and  strengthened  capillary  tube  pressure  vessel- 
sample  containers.  The  shear  viscosity  is  expected  to  show  an  increasingly  well 
resolved  minimum  at  around  150  MPa  as  the  temperature  decreases. 

The  bulk  viscosity,  which  can  be  obtained  from  combining  shear  viscosity 
and  ultrasonic  absorption  measurements,  has  been  given  less  attention  for 
experimental  reasons,  but  its  behavior  down  to  -20°C  is  known  to  be  similar  in 


character  to  that  of  the  shear  viscosity. 


(68,70) 


DIFFUSIVITY  The  most  fundamental  of  the  transport  properties  is  the  self 

diffusion  coefficient  since  no  external  stress  is  required  to  manifest,  or 

measure,  the  property.  Data  have  been  obtained  below  0°C  both  by  tracer, 

(42) 

and  NMR  spin-echo  methods,  the  latter  appearing  more  reliable  and  extending 

o 

to  -31  C.  Combined  with  high  pressure  measurements  performed  utilizing  the 
(95) 

latter  method,  the  diffusivity  shows  the  positive  pressure  dependence 

expected  from  the  viscosity  behavior,  though  the  maximum  on  the  high  pressure 
side  is  not  very  clearly  defined.  Self-diffusion  measurements  on  emulsion 
samples  are  possible  in  principle  using  a  pulsed  gradient  NMR  spin-echo  tech¬ 
nique,  though  the  interval  between  pulses  must  be  kept  extremely  short  in  order 
for  free  diffusion  conditions  to  be  maintained  during  the  experiment.  Since 
viscosity,  but  not  self-dif fusivity,  diverges  near  the  common  critical  and  con¬ 
sulate  points,  the  comparison  of  diffusivity  and  viscosity  down  to  the  low 


temperature  limit  in  the  case  of  water  should  be  a  matter  of  some  interest. 

Some  detail  on  the  diffusion  process  is  available  from  incoherent  quasi¬ 
elastic  neutron  scattering  ezpriments  which  were  recently  extended  to  -20°C  by 
Chen  et  al'  This  work  clearly  resolved  two  Lorentzian  components  to  the 

(instrument-resolution-function-corrected)  quasi-elastic  spectrum.  The  broad 

component  characterized  by  a  halfwidth  P.  of  (3.3  ps.)  1  at  294.4  K  was  associ- 

0 

ated  with  rotational  relaxation  (the  NMR  rotational  correlation  time  at  the 
same  temperature  is  3.2  ps.).  The  sharp  (s)  component  was  linked  to  proton 

translational  diffusion  and,  using  a  jump  diffusion  model  for  the  analysis,  a 

fi  2  2 

characteristic  jump  length  was  extracted  from  a  plot  of  j  ^/DQ  vs.  Q  .  The 
mean  jump  distance  was  1.23  ft,  independent  of  temperature  between  253  and 
273  K.  This  corresponds  to  the  proton  moving  from  one  hydrogen  bond  site  to 
another  via  a  partial  molecular  rotation  across  the  tetrahedral  angle. 
According  to  the  model  this  motion  would  occur  on  average  every  1.1  ps  at  298  K 
and  every  5.7  ps  at  -20°C,  the  temperature  dependence  being  that  of  the  diffu¬ 
sion  coefficient. 

This  characteristic  time  and  its  implied  temperature  dependence  is  longer 
than  that  associated  with  the  breaking  hydrogen  bond  on  the  basis  of  depolar¬ 
ized  light-scattering  (discussed  at  the  end  of  this  section)  though  the  two 
processes  would  seem  to  be  closely  linked.  It  may  be  pertinent  to  observe  that 
jumps  of  the  distance  described  only  make  sense  in  the  context  of  a  highly 
structured  system,  and  a  proton  that  'disconnects'  in  a  non-structured  micro¬ 
region  will  behave  in  a  less  regular  manner. 

NMR  RELAXATION  The  deepest  penetration  of  the  supercooled  regime  achieved  for 
any  mass  transport  or  relaxation  property  has  been  for  the  case  of  spin-lattice 
relaxation  times,  for  which  sample  dimension  is  not  a  problem,  and  the  emulsion 
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sampling  method  may  be  used.  Normal  pressure  studies  on  emulsions  were  first 


(48  97) 

made  by  Hindman  et  al .  '  and  the  technique  haa  aince  been  exploited  with 


great  success  by  Lang  and  Ludemann 


(34,98-100) 


who,  using  Tamada's  HF-strengthened  glass  for  x, 


sample  containment,  have  explored  the 


pressure  regime  up  to  300  MPa,  and  pushed  the 

o 


measurements  down  to  as  low  as  -86  C. 

M 

The  results  of  Lang  and  Ludemann’ s 


1  2 

measurements,  conducted  variously  on  H,  H, 


17 


and  0  nuclei,  have  been  most  enlightening. 


Most  spectacular,  perhaps  is  the  observation  <»'  3 


of  a  minimum  for  protons  and  deuterons  un¬ 


der  applied  pressures  of  the  order  of  200  MPa 


near  -80  C.  The  implication  of  the  minimum  is 
that  water  under  these  conditions  is  in  a 
highly  viscous  condition. 

Orientation  correlation  times  for  the 
water  molecule  have  been  deduced  from  such 


measurements  for  all  3  nuclei,  and  their 


concordance  is  demonstrated  in  Figure  6  for  Figure  6.  Test  of  Eq.  (2)  for  T2 

obtained  from  iH(^S),  ZH(W)  ,  and 
0(toO)  nuclear  spin-relaxation  time 


data 


'V' 


the  pressure  range  3-150  MPa,  plotted  in  the 

Eq.  (2)  form.  The  support  given  this  simple  equation  by  these  data,  which 

extend  over  some  3  orders  of  magnitude  in  r^,  is  rather  impressive.  The 

variation  of  T  with  pressure  closely  mirrors  that  of  T.  in  the  same  way  as 
t  lx  * 


found  for  the  compressibility.  The  relationship,  for  both  H^0  and  D^O,  is 


shown  in  Fig.  7. 
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For  pressure  above  150  MPa,  Lang  and  Lndemann  find  Eq.  (2)  leaa  satisfac¬ 
tory,  or  at  least  indistinguishable  in  its  ability  to  describe  the  behavior 
from  the  familiar  Vogel-Tammann-Fulcher  equation 

X  -  A  exp  +  B/(T-T  )  W 

x  o 

which  describes  the  transport  properties  of  most  viscous  liquids.  Thus  the 
growth  of  the  density  and  entropy  fluctuations  which  seem  responsible  for  the 
low  pressure  anomalies  of  water  is  inhibited  by  '  increasing  pressure  and  their 
effects  are  displaced  to  lower  temperatues  where  their  influence  on  the 

transport  properties  is  finally  masked  by  the  cooperative  relaxation  phenomena 

-T- 

characteristic  of  the  viscous  regime.  ,K)  ,  ' 


In  the  high  pressue  regime,  the  physi-  j  j  i 

\  h2o 

cal  significance  of  the  VTF  equation  para-  240j — .  j 

meters  can  be  critically  examined. ^4,41,65) 


The  pre-exponent  has  the  role  of  an  attempt 

200 

frequency  and  might  be  expected  to  correlate 
with  the  vibrational  frequency  observed  in  >80 
the  IK  spectrum,  hence  to  exhibit  an  isotope  Q 
effect  which  is  predictable  viz.. 


j  S'-  v 

1  >S 
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.  pIMFo) 

The  extent  to  which  this  expectation  is  Figure  7.  Comparison  of  the  sin¬ 
gular  temperatures  Ts  for  H2O  and  D2O 

borne  out,  and  the  additional  expectation  -  obtained  by  Eq .  (2)  analysis  of  re¬ 
laxation  time  data,  with  the  homogene- 

that  the  T  parameter  for  D  O  fall  some  5°  ous  nucleation  temperature  Th-  Also 
0  2  shown  are  the  "ideal"  glass  tempera- 

above  that  for  ILO  as  is  found  from  direct  tures  T0  obtained  from  Eq.  (3)  fits 
2  to  data  above  150  MPa.  (After  Lang 

glass  transition  temperatue  measurements  on  and  Ludemann,  ref.  34,  reproduced  by 

permission) . 

H^O  and  D^O  solutions  -  is  upheld,  may  be  judged  from  the  following  table  of 

*  (34) 

best  fit  parameters,  presented  by  Lang  and  Ludemann. 
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Tjhk  I  Parameter*  of  ihr  VTh -equation  oht.imcU  from  the  least  squares  i tiling  of  r^io  t<\  (■*>  in  •standard  deviation) 


p  jMPa)  r..±:|k|  B±n  |k)  (r.irrjxM)  4  M 

D  O  H  O  “  ‘  ~  •'  '' 


Bin  |kl 
D;U  110 


(r„2m  x  M)  4  js) 

l)  O  H  .O 


652  ±5 

6  6*0  J 

46*0.3 

62.1  *  7 

6  liOS 

4.6*0  3 

5KI  ±5 

7.0 £  0  6 

5  V*  0  3 

C  nrr  -coeff 
0.0  H;0 


The  value  the  Eq.  (4)  Tq  Is  in  accord  with  thermodynamic  estimates,  as 

(41) 

best  they  can  be  made  at  the  moment.  These  results  lend  support  to  the 

notion  discussed  later  that  the  transport  properties  of  water  (and  its  dilute 
solutions)  be  analyzed  as  a  superposition  of  anomalous  components  of  Eq.  (2) 
form  and  a  background  component  of  Eq.  (4)  form. 


OTHER  RELAXATION  AND  DYNAMICS  STUDIES  Information  on  the  reorientational 
motions  of  the  water  molecules  can  also  be  obtained  from  the  dielectric  relaxa¬ 
tion  phenomenon  which  has  been  studied  by  various  workers^®^’*^  using  dif¬ 
ferent  techniques  penetrating  as  far  as  -21°  into  the  supercooled  regime.  The 

most  recent  of  these  studies  and  also  the  one  to  achieve  the  lowest  tempera- 

a 

tures  is  that  of  Bertolini  e_t  al^*^  who  used^coaxial  resonator  technique 
requiring  samples  of  only  0.1  cc,  and  achieved  good  accuracy  in  their  data. 
These  workers  were  able  to  show  that  the  single  relaxation  time  characteristic 
of  water  in  its  normal  regime^^  is  retained  down  to  the  lowest  temperatures 
studied.  This  phenomenon  is  somewhat  uncommon  amongst  molecular  liquids  and 
seems  restricted  in  general  to  the  monohydric  alcohols.  Its  occurrence  in  an 
associated  multi-laterally  hydrogen  bonded  liquid  like  water  is  unexpected  and 
suggest ive . 

The  dielectric  relaxation  time  and  the  structural  relaxation  time  deter¬ 
mined  by  frequency-dependent  mechanical  relaxation  studies  (e.g.,  ultrasonics, 
Brillouin  scattering)  are  usually  found  to  lie  close  to  one  another  in  the  case 
of  associated  liquids. If  this  were  the  case  for  water  then  the  Brillouin 
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scattering  measurements  conducted  by  Teizeira  and  colleagues^^'^^  would  be 
expected  to  show  a  broadening  of  the  peaks  at  around  -20°C  where  the  dielectric 
relaxation  time  is  found  from  Bertolini  <^t  al's  work,  to  be  *»«”*,  the  inverse 
of  the  frequency  of  the  90°  excited  (or  annihilated)  phonon  frequency.  How¬ 
ever,  there  is  no  indication  of  such  broadening  in  the  measurements  reported  to 
date,  suggesting  that  the  longitudinal  relaxation  time  is  at  least  an  order  of 
magnitude  shorter  than  the  dielectric  relaxation  time,  as  estimated  by  Slie  et 
al.dO^)  An  explanation  for  this  difference  is  not  yet  at  hand.  Presumably  it 
would  be  of  considerable  help  in  understanding  the  microscopic  fluctuations 
which  all  the  other  measurements  indicate  are  present  in  this  liquid.  There  is 


particular  interest  in  the 
possibility  of  two  relaxation 


T/°C  [HjO  only] 
i2 _ 20 _ 0 _ -is 


time  scales  in  the  low  tem¬ 


perature  liquid,  since  the 


020 

l/T.  (2H)/sec'20- 


time  scale  for  critical  flue-  Tj/cP 

tuitions  in  the  usual  critical 
phenomena  are  known  to  be 
longer  than  those  characteris-  H2Q 


tic  of  the  'normal'  liquid. 


Even  the  dielectric  relaxation 


data  at  0  C  give  some  sugges¬ 
tion  of  a  second,  very  weak 


relaxation  at  a  frequency 


l0V/seccm2  W'J*" 

or  05  ^  " 


10"  T0/sec  2°|" 


10'%/sec 


lO/TrfWsec1  ’  /f 

or  0  2  ■/ 
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about  10  times  higher  than  the  1 0/T.( 1  d/sec '  I 


primary  relaxation  (R.  H. 


log!  T/  228- 


Cole,  private  communication). 

Figure  8.  Tests  of  Eq .  (2)  for  various  trans- 
The  conformity  of  the  pert  properties  and  relaxation  tines  in  super¬ 
cooled  water  and  D^O  ;  from  ref.  reproduce.: 

various  transport  coefficients  with  permission j . 


and  relaxation  times  measured  for  H^O  and  D^O  at  normal  pressure  to  Eq.  (2)  is 
demonstrated  in  Figure  8,  the  singular  temperature  favored  by  the  data  being 
228  K. 

The  shortest  time  scale  for  a  dynamic  process  in  liquid  water  is  that 
implied  by  a  second  broad  Lorentzian  observed  in  the  spectrum  of  polarized 
light  scattered  by  non-propagating  refractive  index  fluctuations  in  water. 
This  Lorentzian  exhibits  a  half-width  which  at  0°C  corresponds  to  a  relaxation 
time  for  the  index  fluctuation  of  0.3  ps.  It  was  first  described  by  Montrose 
et  al^10^  in  1974  and  was  recently  studied  down  to  -17°C  by  Conde  and 
Teixeira^*^  using  a  capillary  cell. 

As  in  the  quasi-elastic  neutron  scattering  study  discussed  earlier,  the 
second  Lorentzian  components,  in  this  case  the  narrow  one,  can  be  associated 
with  a  relatively  long  time  constant  process  associated  with  rotational  diffu¬ 
sion;  its  time  constant  and  non-Arrhenius  temperature  dependence  correlate  well 
with  those  of  dielectric and  NMR^^  100)  studies.  The  short  time  process 
on  the  other  hand  has  a  much  smaller  temperature  dependence,  which  as  best  can 
be  discerned  from  the  data  to  hand,  follows  an  Arrhenius  law.  Furthermore  its 
pre-exponent  has  a  value  of  5.6  x  10  *^sec  corresponding  (2nrc)  *  to  the  broad 
IR  librational  band  frequency  of  ~700  cm  *  ,  consistent  with  a  process  involv¬ 
ing  frequent  attempts  to  cross  a  librational  barrier  with  occasional  successes 
in  which  tho  proton  'flicks'  through  the  tetrahedral  angle  to  form  a  'new' 
H-bond  at  the  expense  of  an  'old'  one.  This  elementary  restructuring  process 
is  considered  to  proceed  incessantly  and  without  need  of  cooperative  involve¬ 
ment  of  other  molecular  motions  (in  fact  a  temporarily  static  environment  is 
implicit  in  the  model)  -  hence  the  Arrhenius  behavior.  It  is  quite  likely  that 
the  same  local  'switching'  motion  is  executed  and  reversed  several  times  before 
a  structural  rearrangement  occurs.  The  diffusion  of  the  defects  in  the  local 
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quasi-lattice  which  are  created/eliminated  is  presumably  the  process  which 
allows  rotational  diffusion  and  viscous  flow  to  occur.  The  time  constant 
obtained  fom  this  measurement  is  regarded  as  the  'hydrogen  bond  life¬ 
time.  '  ^*05,106)  n£)t  ciear>  however,  why  this  motion  is  not  seen  as  a 

separate  contribution  in  NMR  studies. 

SOLUTION  STUDIES  If  instead  of  studying  pure  water,  we  examine  the  properties 

of  dilute  solutions  in  which  the  foreign  species  can  serve  as  a  probe  for  the 

dynamics  of  the  solvent  structure  in  which  they  move,  we  can  in  principle  gain 

a  considerable  increase  in  the  precision  with  which  the  transport  properties 

may  be  characterized.  For  instance,  Speedy^1®^  has  recently  obtained  data 

precise  to  1%  on  the  electrical  conductivity  of  solutions  of  EC1  and  HC1  in 

water  using  fine  capillary  conductance  cells  which  permit  an  extension  of  the 

data  to  -30°C  at  pressures  up  to  50  MPa.  This  study  considerably  extends  an 

earlier  penetration  of  the  supercooled  regime  by  Takisawa  et  al.^^**)  Speedy 

found  that  his  low  temperature  data  at  the  lowest  concentrations  could  be 

fitted  to  Eq.  (1)  with  the  same  parameters  which  were  yielded  very  precisely  by 

a  fit  of  the  more  precise  and  accurate  literature  data  for  temperatures  above 

0°C.  The  best  fit  of  T  from  these  studies  is  221.*,  E  at  the  lowest 

s 

concentrations.  The  consistency  of  these  results  is  an  important  indication 
that  the  supercooled  water  data  obtained  in  fine  capillaries  is  in  no  way 
influenced  by  contamination  effects,  surface  effects,  etc.  provided  excessively 
fine  capillaries  are  avoided. 

An  interesting  feature  of  the  conductivity  fits  is  that  the  exponent  best 
fitting  the  data  is  significantly  lower  than  those  characterizing  the  scopes  of 
the  plots  in  Figure  8.  In  the  case  of  dilute  HC1  solutions.  Speedy  et  al/109* 
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have  osed  the  chloride  ion  conductivity  derived  from  the  EC1  data  to  extract  a 

conductivity  vs.  temperature  relation  for  the  proton  X.  .  The  proton,  of 

coarse,  is  restricted  to  a  peculiar  mode  of  transport  which  links  its  behavior 

closely  to  the  hydrogen  bonded  network  of  the  solvent.  In  this  case,  the  Eq. 

(2)  with  an  exponent  of  anity  describes  all  the  data,  and  yields  perhaps  more 

convincingly  than  in  any  previous  study,  the  value  of  the  singular  temperature, 

T  =  227.  4E. 
s 

(91) 

In  related  work.  Speedy  et  al  have  tested  the  effects  on  the  dilute 

KC1  conductance  of  ethanol  additions  while  exploring,  for  transport  properties, 
the  solute  effect  described  earlier  in  this  review  in  connection  with  thermo¬ 
dynamic  properties. The  interesting  result,  found  also  by  Halfpap  and 
( 89) 

Sorensen,  see  below,  emerges  that  although  the  density  fluctuations,  judged 

from  the  behavior  of  the  compressibility,  are  rapidly  damped  out  by  such  addi- 
(71) 

tions,  the  anomalies  in  the  conductivity  are  enhanced.  There  is  slight 

evidence  for  an  actual  increase  in  the  singular  temperature.  Certainly  the 
latter  is  not  depressed  as  it  is  by  additions  of  ionic  salts. In  analyz¬ 
ing  these  data  Speedy  et  al  introduced  the  notion  of  background  terms  in  the 
transport  equation.  Choosing  a  VTF  equation  (Eq.  (4)  above)  form  for  the 
latter,  they  obtained  data  fits  to  the  expression 

l/o  -  (T/Ts-1)A  exp[-B/(T-160)]  ^ 

in  which  the  exponent  in  the  singular  term  of  Eq.  (2)  now  has  the  value  of 
unity.  The  implication  is  that  exponents  differing  from  unity  in  fits  to 
Eq.  (2)  reflect  the  presence  of  background  terms,  as  was  the  case  for  thermo¬ 
dynamic  properties. ^  This  would  be  consistent  with  the  notion^*1^  that  both 
transport  and  thermodynamic  properties  are  analytical  functions  of  temperature 
about  the  singular  temperature,  an  assumption  which  yields,  by  a  Taylor  expan¬ 
sion,  a  leading  term  in  (T/Ts”l)  ,  i.e.  7  ■  1.0. 
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A  justification  for  this  form  of  background  is  provided  by  an  analysis  of 

binary  solution  data  on  the  system  H^O  +  hydrazine  performed  by  Oguni.^^^ 

Following  the  solution  extrapolation  principles  developed  by  Oguni  and 
(87) 

Angell,  viscosities  of  hydrazine  +  water  solutions  measured  by  capillary 

viscometry  at  temperatures  ranging  down  to  -35°C  were  extrapolated  to  zero 

water  content  to  obtain  a  viscosity-temperature  relation  for  'normal'  water. 

The  'normal'  viscosity  vs.  temperature  relation  for  water  was  found  to  conform 

to  the  VTF  equation  with  best  fit  parameters  similar  to  those  found  from  the 

high  pressure  studies  relaxation  time  in  Table  1  above,  viz.,  T  *  143  E, 

o 

B  =»  409  E  .  The  fit  using  the  parameters  extrapolated  from  the  high  pressure 

studies,  Tq  *  128  E  ,  B  *  520  E  is  only  slightly  poorer,  and  the  latter  value 

is  more  consistent  with  an  independent  estimate  based  on  a  thermodynamic 

(41) 

vanishing  entropy  analysis. 

The  use  of  solvent  effects  for  diagnosing  the  structural  characteristics 

associated  with  anomalies  in  pure  water  has  also  been  introduced  by  Halfpap  and 
(89) 

Sorenson,  who  used  both  ordinary  capillary  viscometry,  and  light-scattering 
from  latex  spheres,  to  observe  the  effect, of  ethanol  on  the  one  hand  and  hydra¬ 
zine  on  the  other,  on  the  viscosity  of  water  at  low  temperatures.  Here  the 
interesting  observation  was  that  the  completely  hydrophilic  molecule  hydrazine 
reduced  the  viscosity  of  water  quite  sharply,  while  ethanol  with  its  hydro- 
phobic  ethyl  group  strongly  enhanced  the  viscosity  (as  it  did  the  resistivity) . 
Thus,  although  hydrogen  bonding  per  se  is  not  destroyed  by  the  addition  of 
hydrazine,  the  special  hydrogen  bonding  patterns  which  are  associated  with  the 
anomalies  in  water  are  dismantled.  By  contrast,  they  are,  if  anything,  en¬ 
hanced  by  the  presence  of  the  hydrophobic  group.  This,  like  the  results  of 
Speedy  et  al,^^  strongly  hints  that  the  cage-like  water  molecule  arrangements 
known  to  form  around  hydrophobic  groups,  are  the  important  structural  elements 
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associated  with  the  anomalous  character  of  low  temperature  water.  This  is 

( 85 ) 

essentially  the  picture  anticipated  in  the  analysis  of  Stillinger. 

Finally  in  this  section  we  mention  some  very  recent  results  by  Easteal  and 
Woolf,  which  although  not  strictly  dealing  with  supercooled  conditions, 

are  so  diagnostic  for  low  temperature  water  structure  that  they  must  be  in¬ 
cluded  in  this  discussion.  Easteal  and  Woolf  introduced  the  pressure  variable 
into  the  study  of  tracer  diffusion  of  guest  molecules  and  found  striking  dif¬ 
ferences  in  the  pressure  coefficients.  Thus,  where  the  self-diffusion  coeffi¬ 
cient  of  water  increases  with  pressure,  leading  to  the  expectation  that  guest 
molecules  would  exibit  a  comparable  dependence,  Easteal  and  Woolf  find  that  the 
diffusivity  of  the  non-hydrogen  bonding  molecule  acetonitrile,  CH^CN,  actually 
decreases  with  increasing  pressure,  a  direct  violation  of  the  Stoke s-Einste in 
equation.  Ethanol  diffusivity,  on  the  other  hand,  showed  an  even  more  rapid 
increase  with  pressure  than  does  water  itself.  The  only  previous  result  com¬ 
parable  to  the  observation  for  acetonitrile  is  that  reported  for  interstitial 
ions  in  open  network  liquid  silicates  where  both  computer  simulation  stud- 
ies^113^  and  experimental  observatiors  on  the  responses  of  calcium  diffusion 
and  viscosity^^  to  pressure  increases  in  closely  similar  liquid  silicates, 
show  the  interstitial  ions  decreasing  in  mobility  while  the  liquid  fluidity  in¬ 
creases.  The  relationship  between  liquid  silica  (and  open  framework  silicates) 
and  water  has  recently  been  emphasized, the  Easteal-Woolf  observations  are 
presumably  to  be  interpreted  as  for  the  liquid  silicate  interstitials.  In 
other  words^  their  observations  suggest  the  existence  in  liquid  wate-  at  low 
temperatures  of  cage-like  structures  through  which  non-hydrogen  bonded  mole¬ 
cules  can  slip  at  low  pressures  but  in  which  the  channels  become  constricted 

and  inhibiting  as  pressure  increases.  The  greater  increase  of  D(ethanol),  on 

(89  90) 

the  other  hand,  implies  as  Sorensen  suggested  '  the  location  of  EtOH 


J 
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molecules  in  the  most  open,  highly  structured  aicroregions  of  the  liquid. 
These  are  the  regions  that  respond  most  rapidly,  according  to  LeChatel ier ' s 
principle,  to  increases  in  pressure,  thereby  releasing  the  guest  molecules  for 
more  rapid  diffusivity  increases,  on  average,  than  for  the  host  molecules. 

Simulations  and  Theoretical  Interpretations 

Before  discussing  directly  and  in  detail  the  various  schemes  within  which 
the  above  described  properties  may  be  interpreted,  it  is  appropriate  to  intro¬ 
duce  some  of  the  additional  general  observations  made  possible  by  molecular 
dynamics  calculations.  Although  further  reference  to  analyses  of  the  output  of 
these  calculations  in  support  of  the  various  theoretical  schemes  will  be  made 
below,  certain  items  of  interest  in  the  sense  of  additional  data  input  should 
first  be  described.  For  the  case  of  the  deep  supercooled  regime,  simulations 
are  not  extensive:  due  to  the  long  structural  equilibration  times  involved,  the 
region  is  out  of  reach  of  all  but  the  most  powerfully  funded  simulationists. 
(To  experimentalists  who  find  themselves  sometimes  irked  by  the  apparent  ease 
with  which  current  simulation  programs  can  reproduce  and  then  interpret  their 
most  painstakingly  achieved  results,  this  situation  comes  as  a  mild  relief.) 

Apart  from  one  of  the  early  runs  in  the  classic  study  of  Rahman  and 

Stillinger^®’**^  -  a  run  which  settled  down,  at  constant  energy,  to  a  final 

o 

temperature  of  -8.2  C  ~  the  only  significant  penetration  of  the  supercooled 
regime  obtained  using  an  accurate  pair  potential  has  been  that  from  NPT  en¬ 
semble  Monte  Carlo  calculations  using  the  TIPS  2  potential  recently  reported  by 
Jorgensen.  The  results  for  25°C  and  75°C  were  in  excellent  accord  with 

experiment.  At  -30°C,  the  intensive  thermodynamic  properties  Cp,  a,  and 
proved  very  slow  to  converge  and,  while  qualitatively  maintaining  a  correct 
trend  and  exhibiting  an  expansivity  zero  (i.e,,  density  maximum),  failed  to 
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register  the  large  increases  observed  experimentally,  presumably  due  to  the  ex¬ 
cessive  equilibration  time  requirements.  Monte  Carlo  calculations  suffer  the 
same  ergodicity  problems  as  molecular  dynamics  calculations,  and  since  the 
relaxation  times  in  experiment  increase  by  more  than  an  order  of  magnitude  be¬ 
tween  25°C  and  -30°C  the  computational  difficulties  can  be  appreciated.  How¬ 
ever,  even  without  complete  equilibration  being  assured,  the  results  at  -30°C 


were  sufficient  to  underscore  the  importance  of  certain  of  the  key  trends  as 
the  supercooled  regime  is  entered.  To  this  author,  the  most  striking  feature 
is  the  rapid  narrowing,  and  peak  value  increase,  in  the  0-H...0  angle,  coupled 
with  the  large  increase  in  4-bondedness  -  a  requirement  for  anomalous  density 
decreases  at  low  temperatures.  These  rsults  indicate  a  way  of  distinguishing 
tetrahedral  liquid  networks  with  'sharp'  liquid  state  density  maxima,  from 
those  with  weak  or  absent  maxima ^H2°'  1^5°  (sharp)  SiO^,  140°, 

(weak)  '-GeOp  133°,  (absent)  . 

At  a  more  qualitative,  but  still  valuable,  level  the  properties  of  a  wide 
region  of  PVT  space  for  water  have  been  mapped  out  in  a  study  by  Tamamoto  and 
co-workers, who  have  paid  considerable  attention  to  the  supercooled 
regime.  These  workers  used  a  particularly  simple  potential,  which  permitted 
very  fast  Monte  Carlo  calculations  without  sacrificing,  it  would  seem,  the 
essential  requirements  for  reproducing  'water-like*  behavior.  The  equation  of 
state  obtained  from  fits  to  state  points  evaluated  over  wide  ranges  of  tempera¬ 


ture  and  pressure,  reproduces  the  experimental  ob«ervations  of  maxima  at 

low  pressures  and  constant  pressure  divergences  in  the  _  and  a  with  T  values 

1  s 


remarkably  close  to  those  yielded  by  experiment. 


(34,40,107) 


A  refinement  of  the  use  of  simulation  which  is  being  developed  to  clarify 


the  meaning  of  'structure'  in  water  (and  other  liquids)  is  currently  being 
described  by  Stillinger  and  Weber.  The  principle  here  is  to  use  the 
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ability  of  computer  simulation  to  quench  instantaneously  a  typical  configura¬ 
tion  of  a  liquid  such  that  it  drops  into  the  nearest  minimum  on  the  3N+1 
dimensional  potential  energy  surface  conceptualized  by  Goldstein  for  his 
extended  discussions^^)  0£  viscous  liquids,  the  glassy  state  and  the  Kauzmann 
paradox^^D  some  time  ago.  The  liquid  system,  especially  at  lower 
temperatures  and  longer  relaxation  times  is  conceived  as  spending  most  of  its 
time  vibrating  anharmonically  in  the  vicinity  of  such  minima,  the  depth  of  the 
minimum  depending  on  the  temperature.  Stillinger  and  Weber  have  succeeded  in 
reducing  this  picture  to  a  formal  representation  of  the  partition  function  in 
which  the  vibrational  and  configurational  contributions  to  the  configuration 
integral  are  uniquely  separated.  The  idea  behind  the  Stillinger-Weber 
configuration  quench  simulations  is  that,  in  this  way  the  energies  of  these 
structures,  and  their  characteristic  geometries  can  be  determined  .in  the 
absence  of  vibrational  smearing .  by  examining  such  quenched  states  and  thus  the 
broad  features  of  the  potential  surface  can  be  mapped  and  the  configurational 
part  of  the  partition  function  properly  evaluated.  The  structures  involved  in 
the  cooperative  region  below  — 20°C,  and  the  origin  of  the  cooperative  effect, 
should  become  much  clearer  from  the  examination  of  such  'hidden 
structures '  than  it  can  be  from  the  instantaneous  liquid  configurations  or 
even  vibrationally  averaged  configurations  on  which  current  analyses  are  being 
conducted. 

We  should  note,  in  connection  with  the  validity  of  the  Stillinger-Weber 
endeavor,  that  effectively  instantaneous  quenches  are  routinely  performed  by 
glass  scientists  when  they  anneal  glasses  to  constant  'fictive'  temperature 
over  periods  of  days,  and  then  cool  to  room  temperature.  Relative  to  the 
relaxation  time,  this  cooling  changes  the  temperature  so  rapidly  that  no 


structure  change  other  than  by  vibrational  contraction  can  occur.  This  is 
possible  because  the  system  near  its  glass  transition  temperature  is  already 
confined  to  exploring  such  deep  minima  that  it  has  no  possibility  of  escaping, 
during  the  final  cooling  to  room  temperature,  the  minimum  above  which  it  is 
located  at  the  time  of  quenching.  (At  the  Kauzmann  temperature  it  would,  in 
principle,  be  in  the  single  lowest  minimum  on  the  surface  -  which  it  would  in 
principle  take  an  infinite  time  to  find!)  By  instantaneous  quenches  in 
molecular  dynamics  simulations,  systems  which  exist  in  palpably  fluid  states 
can  be  trapped  in  one  of  (the  nearest  one)  the  much  shallower  and  more  frequent 
minima  characteristic  of  the  higher  levels  of  the  potential  surface.  Such 
quenches  have  been  carried  out  for  hard  spheres  and  liquid  KC1  to  establish  the 
'loose  random  packing'  limit and  they  are  much  needed  for  the  case  of 
water,  particularly  from  supercooled  states  to  help  elucidate  the  'significant 
structures'  problem. 

In  the  meantime  we  examine  the  various  models  which  have  been  proposed  to 
help  understand  the  supercooled  regime,  recognizing  that  a  model  adequate  to 
deal  with  this  extreme  regime  will  probably  have  much  to  say  about  the  stable 
state  of  water  in  which  the  greatest  interest  naturally  resides. 

For  reasons  of  space,  our  treatment  cannot  be  comprehensive,  and  we  will 
emphasize  the  approaches  we  consider  most  profitable  and  interesting. 

At  this  time  there  appear  to  be  two  distinct  classes  of  model  being 
developed: 

(1)  those  which  focus  on  the  statistical  consequences  and  correlations, 
which  develop  from  random  bonding  patterns  subject  to  the  geometrical  restric¬ 
tion  of  a  maximum  of  4  bonds  per  molecule, 

(2)  those  which  seek  to  identify  some  cooperative  principle  which  directs 
the  course  of  bond  pattern  development. 
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Additionally  there  is  some  support  for  a  'heterophase  fluctuation'  inter¬ 
pretation  of  the  phenomena  in  which  increases  of  heat  capacity,  for  instance, 

are  associated  with  production  within  the  liquid  of  ultra-micro  prenucleation 
(9  72) 

crystallites  '  with  approximately  the  enthalpy  per  mole  of  the  solid  phase 

(ice  I).  However,  this  interpretation  will  not  be  discussed  further  here 

because  the  present  author  is  unable  to  understand  why,  if  this  approach  were 

correct,  the  phenomena  approaching  the  homogeneous  nucleation  temperature  are 

not  manifested  when  other  liquids,  or  even  an  80%  water  +  20%  H  O  solution, 

*  * 

are  supercooled. 

Ve  should  note  also  that  the  'background'  liquid  on  which  the  anomalies 

superimpose,  seems  to  be  very  well  described  by  the  Sceat-Rice  random  network 
(23) 

model.  Thescjauthors  suggest  that  the  absence  of  anomalies  in  the  tempera¬ 

ture  regime  -45  to  25°C  is  due  to  the  absence  of  ring  closures  in  their  Cayley 
tree-based  model.  This  excludes  the  possibility  of  large  density  fluctuations 
since  without  ring  closures  the  bulky  structures  characteristic  of  ice  I  and 
the  clathrates  cannot  arise,  and  the  density  maximum  will  not  appear. 

Predominant  among  models  of  Class  1  is  the  4-bonded  molecule  percolation 

( 81 ) 

model  of  Stanley  and  Teixeira  (ST)  which  has  been  much  discussed  recent¬ 
ly.  This  model  examines  the  statistical  consequences  of  increasing 

'bond*  connections  -  defined  by  an  energy  cut-off  criterion  -  which  form  with 
decreasing  thermal  energy,  and  focusses  attention  on  the  adjacency  of  molecules 

which  are  connected  to  4  other  molecules.  Using  an  ice  I,  lattice  and  random 

n 

placement  of  bonds,  ST  show  how  clusters  of  4-bonded  molecules  naturally  form 
and  rapidly  proliferate  as  the  total  fraction  of  possible  bonded  sites  ap¬ 
proaches  unity.  If,  as  seems  reasonable  (and  as  anyway  can  be  established  by 

analysis  of  an  ST2  simulation  configuration  even  without  determination  of 
(122) 

'hidden*  structure  ),  the  volume  per  molecule  in  a  4-bonded  cluster  is 


greater  than  for  molecules  outside,  or  on  the  edge  of,  such  a  cluster,  then  the 
features  of  water  in  the  supercooled  region  can  be  qualitatively  and  (except 
for  the  lowest  temperatures)  almost  quantitatively  accounted  for.  The  trans¬ 
port  property  temperature  dependences  in  this  model  are  assumed  to  be  deter¬ 
mined  by  the  changing  concentration  of  1-  and  2-bonded  molecules,  an  approach 
which  seems  to  this  reviewer  to  be  an  somewhat  arbitrary  and  simplistic 
treatment  of  what  are  surely  complex  and  cooperative  processes,  but  which  has 
the  advantage  of  describing  the  data  with  a  minimum  of  input. 

(123) 

In  a  current  study  of  ST2  water  characteristics  Geiger  and  Stanley 
have  found  that  the  percolation  exponents  for  the  disordered  continuous  system 
are  essentially  the  same  as  those  known  for  random  bond  percolation  on  3-dimen¬ 
sional  ordered  lattice.  This  tends  to  vindicate  the  original  ice  lattice-based 
formulation  of  the  model. 

The  Stanley/Teixeira  model  repreaents  a  controlled  return  to  the  old  and 
discredited  cluster  models  of  water  with  the  important  difference  that  the 
nature  of  the  matrix  in  which  the  'patches'  exist  is  clearly  described  and  is  a 
continuously  and  variously  bonded  assemblage  (since,  as  argued  some  time  ago  by 
Gibbs  et  al^*^^  and  confirmed  by  MD  configuration  analysis,  the  percola¬ 

tion  threshold  for  single-bonded  connectedness  is  well  above  the  boiling 

point.)  Furthermore,  the  bulky  patches  of  the  ST  model  do  exist  in  the  MD  con- 

(122) 

figurations  of  ST2  water  according  to  a  recent  analysis  of  Geiger  et  al . 

and  have  the  dimensions  expected  from  the  correlation  length  of  density  fluc- 

(54) 

tuations  measured  by  Bosio  et  £l's  small  angle  X-ray  scattering  studies. 

It  is  a  feature  of  the  ST  model  that,  while  a  percolation  threshhold  for 
4-bonded  molecules  appears  to  be  imminent  in  the  vicinity  of  -45°C,  no  singu¬ 
larity  in  the  thermodynamic  properties,  such  as  has  been  indicated  by  the  vari¬ 
ous  data  fits  reviewed  in  early  sections  of  this  review,  is  predicted  by  the 
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model.  Rather,  a  maximum  in  the  isothermal  compressibility  is  predicted  to  lie 
at  about  -60°C.  This  is  in  the  vicinity  of  the  temperature  at  which  the  growth 
of  the  'patches'  to  a  dimension  where  nucleation  of  ice  is  probable  (although 
the  4-bonding  statistics  could  be  equally  well  worked  out  on  a  random  network 
or  a  clathrate-hydrate  lattice  probably  with  little  difference  in  outcome.)  A 
reason  for  some  concern  with  the  ST  treatment  is  that  the  distribution  of  1-, 
2-,  3-,  and  4-bonded  molecules  predicted  by  the  model,  while  according  quanti¬ 
tatively  with  the  output  of  the  Rahman  Stillinger  ST2  computations,  seems  to 
accord  equally  well  with  the  distributions  analyzed  from  other  simulated 
waters,  even  Tamamoto  water, though  the  calculated  physical  properties  of 
the  simulated  waters  are  markedly  different. 

To  obtain  quantitative  results } the  ST  model  assigns  an  excess  volume  per 
mole  to  molecules  involved  in  4-bonded  patches,  and  then  uses  the  quantitative 
assessments  of  the  fraction  of  the  molecules  in  such  patches  in  relation  to 
temperature  to  predict  the  temperature  dependence  of  properties  such  as  the 
compressibility.  This  leads  to  the  prediction  of  a  maximum  in  at  deep 
supercooling.  An  alternative  way  of  producing  such  a  result  with  a  single 
parameter  within  a  random  bonding  model  and  with  a  continuous  relation  between 

bonds  per  molecule  and  volume,  is  via  a  volume-cooperative  bond-lattice 

( 82 ) 

description  in  which  the  AV  of  the  simple  two  state  equation  is  replaced 

with  a  AV  /N  where  N  is  the  fraction  of  possible  bonds  broken  at  the 
o  x  x 

temperature  under  consideration  and  AVq  is  the  volume  increment  per  bond  at 

high  temperatures.  This  approach  was  considered  briefly  in  ref.  40,  but  not 

developed  because  it  did  not  account  for  the  singularity  apparently  implied  by 

the  various  water  anomalies  reviewed.  Better  correlation  with  spectroscopic 
(2  49) 

data  '  would  be  obtained  if  the  volume  increment  per  bond  were  to  change 
sign  at  intermediate  extents  of  bonding,  e.g.,  by  writing 

AV  -  -AV  /N  +  2AV  (6 

OX  o 
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The  problem  of  describing  the  existence  of  supercooling  within  the  lattice 

models  is  briefly  considered  by  Neijer  et  al^^^  bat  no  attempt  is  made  to 

describe  the  properties  of  the  supercooled  state. 

Turning  to  Class  2  ^cooperative)  models  within  which  singularities  might  be 

(85  127) 

sought  we  encounter  two  discussions  by  Stillinger.  '  In  ref.  127 


Stillinger  considers  the  effect  of  system  size  on  such  properties  as  melting 
point  and  density  maximum  and  develops  the  partition  function  for  an  N  molecule 
system  in  terms  of  two  order  parameters  one  of  which. V,  describes  the  'ice- 
likeness*  of  the  system  ( later simpl if ied  for  discussion  in  terms  of  the  concen¬ 
tration  of  (bulky)  bicyclic  octamera)  and  the  other  (g)  the  potential  energy 
(+pV)  for  the  molecules.  He  then  constructs  a Y-  0  surface,  and  shows  how  its 
form  implies  the  possibility  of  a  'catastrophe'  dur.&'t  supercooling  of  the 
liquid  with  diverging  fluctuations  in  the  paraae  ~r,  Stillinger  discusses  the 
conditions  under  which  a  liquid  (water)  could  follow  this  path.  The  essential 
feature  is  the  existence  of  an  attraction  between  bulky  (high  ¥  )  elements  of 
the  structure. 

In  ref.  85  Stillinger  discusses  a  variety  of  possible  molecular  configura¬ 
tions  which  could  be  responsible  in  cold  water,  for  locally  bulky  microregions. 
He  notes  that  typically  these  alternative  structures,  which  consist  of  frame¬ 
works  of  'good*  or  unstrained  hydrogen  bonds,  will  have  faces  which  can  match 
to  one  another  also  without  strain,  hence  should  experience,  in  otherwise 
strained  surroundings,  a  net  attraction  for  one  another.  Since  they  are  them¬ 
selves  internally  low  in  energy  they  will  tend  to  both  increase  in  number  and 
to  cluster  as  temperature  decreases.  Thus  the  ingredients  for  an  accelerating 
cluster  growth,  with  an  eventual  divergence  in  cluster  size  at  a  singular  tem¬ 


perature  T(  are  present.  The  formation  of  the  infinite  ramified  cluster  at  the 
singular  temperature  corresponds  well  with  the  notion  of  'long  range  ordering 
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(18) 

of  strong  hydrogen  bonds'  introduced  by  the  author  and  colleagues  to 
rationalize  the  existence  of  a  lambda-like  anomaly  they  felt  was  implied  by 
their  anomalous  heat  capacity  and  proton  chemical  shift  studies,  combined  with 
the  existence  of  vitreous  ice. 

The  question  to  be  decided  is  whether  or  not  the  infinite  cooperative 
cluster  of  the  Stillinger  model  (which  in  the  ST  model  is  the  percolation 
threshold  for  randomly  formed  4-bonded  molecule  clusters  -  the  difference  is 
not  great)  can  be  reached  without  first  encountering  a  limit  of  mechanical 
stability.  This  alternative  possibility  has  been  the  particular  province  of  R. 
J.  Speedy  who  first  introduced  the  idea  in  the  supercooled  water  context  in  a 
joint  paper  with  the  author and  has  aince  developed  it  in  a  particularly 
interesting  and  provocative  paper  in  which  the  properties  of  both  supercooled 
and  superheated  water  are  tied  together  into  a  single  overall  phenomenological 
scheme. We  will  elaborate  on  this  in  a  concluding  section  below,  but 
first  will  continue  the  discussion  of  cooperative  principles  in  water  since  the 
basis  for  a  long  range  correlation  effect  is  necessary  to  the  distinction  be¬ 
tween  water  and  other  liquids  in  Speedy's  overall  scheme. 

For  his  cooperative  effect  Speedy  takes  inspiration  from  Keppler's  1611 

(128) 

exhortation  to  biologists  to  seek  the  origin  of  the  near-ubiquitous  inci¬ 
dence  of  5-sidedness  in  plants  in  the  sap  that  runs  in  their  structures.  In  a 
current  article  Speedy  refines  the  argument  for  geometrical  cooperativity  out¬ 
lined  earlier, and  argues  that  pentagonal  rings  have  a  natural  propensity 
for  self-replication,  and  that  the  structures  they  form  by  their  nature  exclude 
much  volume  to  other  molecular  centers,  i.e.,  tend  to  generate  low  density 
patches  in  the  liquid.  Some  support  for  the  easy  self-replication  notion  may 
be  drawn  from  the  observation  that  gas  hydrates  may  be  formed  with  many  fewer 
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incorporated  guest  molecules  when  freezing  is  rapid,  (D.  Glew,  private  communi¬ 
cation,  1977)  i.e.,  clathrate  structures  in  the  rapidly  quenched  liquid  appear 
to  propagate  more  efficiently  than  does  Ice  1^  (this  leads  to  some  concern 

about  the  conformations  characteristic  of  water  vitrified  in  the  presence  of 

(20  21  24) 

small  hydrophobic  molecule  cryogens^  *  *  ).  In  a  pentagonal  ring  each  pair 

of  molecules  finds  its  components  in  an  eclipsed  conformation  whereas  in  the 
hexagonal  rings  of  Ice  1^  all  pairs  are  staggered.  This  distinction  provides  a 
basis  for  ongoing  (Speedy,  R.  J.  and  Mesei,  K. ,  private  communication)  analysis 
of  computer-generated  configurations  for  trends  in  Ice  I^-like  vs.  pentagonal 
sided  clathrate  structures  during  supercooling  (Ice  1^  has  two  eclipsed  and 
four  staggered  pairs  in  its  hexagonal  ring  elements). 

STABILITY  LIMIT  CONJECTURE  The  basic  idea  here  is  that  the  large  fluctuations 
in  density  and  entropy  encountered  in  the  deep  supercooling  region  reflect  the 
approach  of  the  system  to  a  state  of  mechanical }  as  opposed  to  thermodynamic, 
stability.  At  the  limit,  commonly  known  as  the  spinodal,  the  fluctuations  grow 
without  limit  carrying  the  system  completely  out  of  its  initial  phase  space, 
the  molecules  reassembling  in  some  other  region  which  will  tend  to  be  the  near¬ 
est  stable  phase.  The  best  known  case  is  the  behavior  of  a  superheated  liquid 
at  a  pressure  less  than  the  critical  pressure.  Here  a  small  pure  droplet  of 
water  at  1  atm.  presure,  for  instance,  will  spontaneously  rupture  at  280°C 
because  of  its  close  approach  to  the  spinodal  predicted,  from  the  most  accurate 
equations  of  state,  to  lie  at  324°C  (it  is  actually  an  imprecise  region  of  tem¬ 
perature  smeared  by  fluctuations,  rather  than  a  point,  according  to  current 
( 129)  ( 1 1 0) 

theory  ).  Speedy  has  taken  the  superheated  liquid  spinodal  as  a  point  of 
departure  for  developing  a  global  picture  for  the  mechanical  stability  boundary 
on  water.  In  this  view  the  high  temperature  spinodal  which,  at  **320°C  crosses 
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into  the  negative  pressure  region,  is 
finally  bent  back  towards  positive  pres-  700 

sures  due  to  the  growing  importance  of 


ii-ZOOW^  j 

/  u.«  \ 


the  cooperative  4-bonded  molecule  'or¬ 
dering'  phenomenon. 

By  assuming  that  the  pressure  can 


be  expanded  along  an  isotherm  as  a  Tay¬ 


lor  series  in  P~Pg  about  the  spinodal  at 
Ps,Ps  »  Speedy  obtained  a  simple  equa- 
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tion  of  state.  This  describes  the  pre-  Flour, S  Estima.ee  locus  el  the  limit  ot  stability  lor  water  r,  er  p,. 

snown  In  relation  to  trie  equilibrium  phase  diagram  =  molting  line, 
cisely  known  PVT  properties  Of  water  in  r,  ■  boBng  ine.  CP  *  critical  point)  ano  toe  experiment  kinetic  limits 

Ot  stability  T „  (•)  values  ot  T,  anfl  p,  tor  H,0  from  Tables  I— III 


the  0-100  C,  0.1-100  MPa  region  with 


II 6.  (□!  values  of  7,  -  7  °C  lor  DjO  from  Table  III  (column 
2)  and  Flgixe  4;  (+)  p,(100  °C)  calculated  Irom  ret  20:  (X)  p,( 0  *C) 
and  p,(25  *C)  estimated  by  Yayanos;  (T„)  values  from  Table  I  tor 


surprising  accuracy,  and  yields  P  ,0  as  superheated  M}0  and  Table  III  tor  supercooled  water:  (7„)  Is  ex- 

S  S  trapouted  to  terminate  at  PJ5°C)«* -190  MPa  as  tSscussedxt  section 

2.3.  (  ,  rtf  |  VO. 

best  fitted  parameters.  The  variation  ]>>3  pe.r  rv^x  c  .©  rs  ) 

n  m  x.  .u _ 1-  rri _ „  IlUvt  OV  H  ri  - 


of  Ps  with  T  is  shown  in  Figure  9. 


Ps>  the  limiting  low  density  for  the  liquid  state  turns  out  to  oe  constant. 

in  the  range  0-100°C  at  essentially  the  density  (0.89  p.  ),  of  the  empty 

icein 

Type  I  clathrate  hydrate  structure! 

Tho  extremum  in  P  occurs  at  35°C  (P  =  -212.4  MPa),  which  is  also  the 
s  s 

temperature  of  the  density  maximum  extrapolated  to  that  (negative)  pressure,  as 
thermodynamic  analysis  requires. Similar  behavior  is  predicted  by  extrap¬ 
olations  of  the  N.B.S.  multiparameter  equation  of  state  for  water  to  negative 
pressures  (Lester  Haar,  private  communication:  see  also  ref.  130  for  equation 

of  state  details) .  We  would  point  out  that  water  in  the  vicinity  of  the  P  ex- 

s 

tremum  at  35°C  is  inherontly  schizophrenic,  being  unable  to  decide  whether  its 
density  fluctuations  are  tending  to  take  the  system  in  the  direction  of  col¬ 


lapse  to  the  vapor  (one  large  void)  or  collapse  to  a  clathrate  hydrate  (many 
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small  ordered  voids).  The  behavior  of  vs.  T  in  Figure  27  (P  ■  0.6  kbar)  in 
the  paper  on  simulated  water  by  Yamamoto  and  co-workers, (c.f.  Figure  3  in 
ref.  74)  ~  which  would  collapse  to  two  vertical  infinitesimally  spaced  lines  on 
approach  to  P^,  at  35°  (see  Figure  9  insert)  -  indicates  the  extreme  case  of  a 
thermodynamically  anomalous  fluid  I 
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